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ABSTRACT 

An existing code for calculating axial turbine performance 
using multiple stream surfaces was modified and made to run 
on the equivalent of an HP-1000 computer system. Calculations 
were made for the geometry of a 485 horsepower dual-discharge 
air-drive turbine for both on and off-design conditions. The 
results were compared with available data obtained at off- 
design speeds. Agreement of the flow rate and horsepower to 
within 51 was obtained. 
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I. INTRODUCTION 



A. DESCRIPTION OF THE TRANSONIC COMPRESSOR TEST RIG 

The Transonic Compressor Test Rig at the Turbopropulsion 
Laboratory (TPL) of the Naval Postgraduate School is shown 
schematically in Fig. 1 and consists of the following major 
components : 

1. Air drive turbine. 

2. Air supply system. 

3. Associated piping including throttling valves at the 
turbine and compressor inlets. 

4. Test compressor. 

The drive turbine is a dual-flow axial air turbine with 
50% reaction. The geometry is given in Table 1. The profile 
shapes of the turbine rotor and of the stator blades are 
identical and the blades are of constant section along the 
radius as shown in Fig. 2. The stator has 31 blades while 
the rotor has 32 (to avoid resonant excitation from wake inter- 
ference) . The two parallel stages of the turbine are designed 
for the following output and total inlet conditions: 

Pressure Ratio: 2.8 

Total Inlet Temperature: 640°R 

Flow rate: 10.85 LBM/SEC 
Horsepower: 485 HP 
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The compressor presently under test is a transonic single 
stage, axial flow compressor. It is instrumented for measure- 
ments of torque, mass flow rate, stagnation temperatures and 
pressures, case and hub wall pressures, and for unsteady 
pressure measurements in the flow field and at the walls. 

The Air Supply System incorporates an electric motor- 
driven multi-stage axial flow compressor manufactured by 
Allis -Chalmers . It can presently supply up to 12 Ibs/sec of 
air at 3 atmospheres, at temperatures between 560°R and 660®R. 
The compressor is rated at 1250 HP and has a controlled 
variable speed drive. 

B. STATEMENT OF THE TASK 

The Transonic Compressor Test Rig was designed to pro- 
vide the means for obtaining experimental data in fundamental 
compressor phenomena. Following the present experiments, an 
experiment to investigate the onset of supersonic unstalled 
blade flutter is planned which would involve replacing at 
least the present compressor rotor by a rotating cascade of 
flat-plate blades. Such a rotor would not be able to produce 
the pressure ratios required to pump the required flow rates 
through the system. Therefore, it has been proposed, that 
a turbocharger compressor be fitted in series with the 
rotating cascade to provide the required flow through it. 

The turbocharger would also be driven using air from the 
Allis -Charmers air supply system. 
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In order to evaluate the feasability of the turbocharger 
installation, it is necessary to determine the mass flow rate 
required by the drive turbine to drive the test compressor 
at a given power and speed. The remaining air to drive the 
turbocharger turbine is then known and the selection of a 
commercially available turbocharger suitable for this 
application can be made. 

Thus, the performance of the air drive turbine must be 
known over the complete speed range. Of particular impor- 
tance, are the required mass flow rates for given values of 
horsepower. The problem, therefore, is to obtain the turbine 
performance map for all pressure ratios and speeds . 
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II. APPROACH 



A. BACKGROUND 

A search of the most recent literature revealed a number 
of analytical methods for the calculation of turbine off- 
design performance. The majority of these used in a finite 
element approach but little information on the relative suc- 
cess of these methods in practice was available. Two alternate 
methods, both used at the Turbopropulsion Laboratory, were 
those of M. H. Vavra and E. Macchi. Each was examined in 
detail . 

The method of Vavra, given in Ref. [1] is a one-dimen- 
sional (meanline) approach using mathematical modelling and 
experimental data to express flow angles and losses. It is 
primarily a method to design turbine blading but may also be 
used to predict turbine performance for a given set of gas 
inlet and operating conditions when the blading geometries 
are specified. It is assumed that the axial velocity is 
constant along the blading from hub to tip. Vavra states 
that this assumtion is reasonable for blading in which the 
tip-to-hub ratio is equal to or less than 1.15. The ratio 
is 1.312 and 1.424 for the drive turbine stator and rotor 
blading respectively. It was thought therefore, that the 
method of Macchi might yield more accurate predictions. 
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Macchi's method is given in Ref. [2]. The method, imple- 
mented by Macchi in a computer program written for the IBM 360, 
was an extension of the work done by R. Eckert [Ref. 3] and 
R. Harrison [Ref. 4] . Eckert wrote a program, following a 
simplified three-dimensional analysis, which could be used to 
predict the performance of a single-stage axial flow turbine. 
Harrison improved the program by modifying the analysis to 
take into account streamline curvature. Both programs were 
based on the three-dimensional method developed by Vavra in 
Ref. [5]. Macchi's principle improvements to the program were 
to introduce the choice of various methods to calculate gas 
outlet angles and loss coefficients. Two methods of calcula- 
ting gas outlet angles are included; those of Ainley and 
Mathieson [Ref. 6] and Traupel [Ref. 7] . Five methods for 
calculating the loss coefficients can be selected; those due 
to Ainley and Mathieson [Ref. 6] , Dunham and Came [Ref. 8] , 

Balje [Ref. 9], Lonherr and Carter [Ref. 10] and Traupel 
[Ref. 7]. 

Macchi's computer program, as documented in Ref. ] 2] , was 
selected for performance predictions of the drive turbine. 

It should be noted that no card deck of the program was avail- 
able, and no results of using the program were available 
other than those included in Ref. [2]. 

B. ANALYSIS 

The method requires the following assumptions; 

1. There are an infinite number of blades in each blade 
row so that blades downstream do not affect upstream conditions. 
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2. The flow is axisymmetric at locations where the 
equation of motion is solved. 

3. The flow is steady and adiabatic. Thus, the total 
enthalpy through the stator remains constant along a stream- 
line and the relative total enthalpy through the rotor remains 
constant along a streamline. 

4. All equations are solved at between blade row loca- 
tions. Increases in entropy occur in the blade row upstream 
of the stations where equations are solved and the entropy 
change along a streamline between blade rows is zero. 

5. The boundary layers on the turbine casing are not 
accounted for. 

The method of solution is as follows: 

1. Assume initial radial positions of the streamlines. 

2. Obtain the axial velocity distribution by solving 
the equation of motion at the stator outlet. The velocity 
distribution into the stator is assumed to be axial, and 
uniform 

3. Obtain stator loss coefficients. 

4. Check overall continuity and adjust the inlet Mach 
number as necessary. 

5. Check the between-streamline continuity, and adjust 
streamline radial positions as necessary. 

6. Repeat this process for the rotor. 

7. Re-cycle all the above calculations, accounting for 
streamline curvature, and repeat until convergence is reached. 



18 



C. METHOD OF SOLUTION 

The computer code written by Macchi was originally run 
on the IBM 360 computer. The program consisted of a deck of 
over 2000 program cards plus over 60 data cards. Since the 
deck could not be located, it was necessary to re-type the 
program from the listing in Macchi 's paper. However, since 
the IBM 360 computer was soon to be replaced in the period 
in which the work was to be carried out, an alternate 
computer was sought. 

The HP-1000 series mini-computer located at TPL was 
selected for two reasons. First, the machine used FORTRAN 
as did Macchi.'s program. Secondly, it would be a benefit 
to TPL to have the program immediately available on the 
laboratory computer. 

The first steps were to analyze Macchi 's program, in 
detail, and then to run it using his example input/output. 

In analyzing the program it became obvious that the computer 
program listing given in Ref. 2, was not the one used to 
obtain the listed output. Numerous discrepancies were found 
in the listing, some of which would have prevented the pro- 
gram from running; others would have caused incorrect results 
to be obtained. A listing of these discrepancies is con- 
tained in Appendix E. When the program was understood and 
flow-charted, it was keyed-in at the HP-1000 computer ter- 
minal. However, modifications were required to accomodate 
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the program within the mini-computer disc-based operating 
system. 

D. MODIFICATION TO THE COMPUTER CODE 

Since there was no card reader, variable input data such 
as turbine speed had to be entered using data or specification 
statements . This contributed in part to the most difficult 
problem, that of program size. The HP-1000 mini-computer uses 
a disc with a storage capability of 19.5 mega-bytes. However, 
the machine memory is only 124 K Bytes, of which only 29 K 
Bytes is available to a programmer. Also, the available 
memory is divided up, or partitioned into two 18 K and one 
11 K partitions, so that no single program can exceed 18 K. 

It was estimated that Macchi's program was over 100 K. So 
it was clear that the program would have to be modified if 
it were to run on the mini-computer. 

The first modification was to remove all subroutines 
from the program that were not actually used. It will be 
recalled that Macchi's program contained five methods for 
calculating loss coefficients and two methods for calculating 
gas outlet angles. It was decided that only the Traupel 
method of calculating loss coefficients would be retained. 
Traupel was selected for two reasons. Firstly, it was the 
method used by Macchi in his example calculations and there- 
fore the modifed program should still reproduce Macchi's 
results. Secondly, the method of Traupel is widely respected. 
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The method of calculating gas outlet angles was totally 
changed. Neither Ainley and Mathieson [Ref. 6] nor Traupel 
[Ref. 7] was used. Both methods required prohibitively 
large sections of computer code. The method selected was 
that of Vavra [Ref. 1] . 

Use of Vavra 's method greatly simplified the program 
because this method predicts gas outlet angles independently 
of the inlet Mach number. Macchi's approach was to use 
Traupel 's method which is dependent on the Mach number of 
the flow into the blade . 

The above simplifications reduced the program size from 
2257 lines to less than 1800 lines. However, this was still 
too large and the program could not be loaded without over- 
flowing the memory. 

The solution to the problem was found in program 
segmentation. In this process, the computer code is divided 
into a main program and several segments. Each segment is 
a "piece" of the original program. The segments are indi- 
vidually compiled and loaded. However, the segments are 
placed into memory only as they are needed to execute the 
overall program. Thus, a very large program can be made 
to run in the available 18 K partition. Since the present 
program was not originally intended for a mini - computer , 
segmentation was not straight forward. The method finally 
arrived at is detailed in Appendix C. Basically, the main 
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program consists of all the subroutines, while the three 
other segments contain coding which enables program flow to 
proceed in a logical manner. 

Successfully segmented, the program was run using Macchi's 
input. An output was obtained which agreed almost exactly 
with Macchi's results. All output quantities were within II 
of Macchi's quantities. The differences were, in all proba- 
bility, due to the different method of calculating gas outlet 
angles . 

After verifying Macchi's program, the drive turbine geo- 
metry was input and the program was run for a given set of 
operating conditions. The results are discussed in the follow- 
ing section. Note; The "verification" of Macchi's program 
amounted to verifying that the computer code now loaded into 
the HP- 1000, was indeed Macchi's code. It was not known 
whether Macchi's output data were a good or bad prediction 
of performance since they were not compared with test results. 
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III. RESULTS OF AXIAL TURBINE PREDICTIONS 



A. USING BOTH COMPLETE AND MODIFIED PROGRAMS 

The drive turbine geometry was input and the following 
solution flow path was selected: 

1. Stator and rotor loss coefficients were functions of 
pressure ratio. 

2. The blockage factor, C*, used in the equation of con- 
tinuity was equal to the total loss coefficient. 

Four operating points were selected to test the validity 
of the program. Three were off-design points at which measured 
data were available and the fourth was the design point itself. 
Table II contains details of the selected test points for 
Run 1 . 

The program variables were then changed and the following 
new solution flow path was selected: 

1. Stator and rotor loss coefficients were those calcu- 
lated by Traupel’s method. 

2. The blockage factor, C* > i^as equal to the profile loss 
coefficient . 

After reviewing the results of Runs 1 and 2, a further 
modification was made to the program. The original program 
contained a subroutine which checked between-streamline con- 
tinuity. If the total mass flow rate at the stator and rotor 
exits was not evenly divided between the five streamlines , 
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the radial positions o£ the streamlines were adjusted and all 
steps were recalculated using the new streamline positions. 
Hence, for Run 3, a subroutine was removed and the main pro- 
gram was modified so that between-s treamline continuity was 
not examined. 

B. COMPARISON WITH MEASURED DATA 

The results of Run 1, 2, and 3 are tabulated in Table III. 
Run 1 showed predictions of mass flow rate which departed about 
6 % from the measured data. However, the horsepower predictions 
were off by as much as 16.17%. Furthermore, the computer pro- 
gram was unable to reach a solution for the design point. 

Run 2 produced worse results as is evident from the table. 
Again, the program was unable to converge to a solution at 
the design point. 

Run 3 produced more acceptable data. Additionally, con- 
vergence to a solution was noticeably faster and a solution 
was obtained at the design point. Because of this, the method 
used in Run 3 was used to map the drive turbine performance. 

The computer program used to obtain the results of Run 3 is 
described in detail in Appendix A and is listed in Appendix G. 
The results of Run 3 are shown plotted in Figures 3 through 8. 

To obtain the plots in Figures 3 and 6, a value of the 
total inlet temperature was approximated by the method of 
Vavra as contained in Ref. [14]. It was assumed that the 
static turbine discharge temperature should not be less than 
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45°F (505°R) . This corresponds to the approximate temperature 
at which condensation of moisture in the air, assuming 100% re- 
lative humidity, will occur. The inlet temperature was given by 



Total 


Inlet Temperature = Statle Outlet Temperature 

1 - h [1 -ci ) 

s'- Y 


where p , 
s ’ 


the total-static turbine efficiency was assumed 



to be 81%, and 6.^ , vras the total to static pressure ratio. The 
total inlet temperature corresponding to each pressure ratio 
is given in Table IV. 

The computer output corresponding to each point on Figures 
3 through 8 is contained in Appendix F. Only one side of the 
dual flow turbine was analyzed, thus, the resulting printed 
values of horsepower, referred horsepower, moment, referred 
moment, flow rate and referred flow rate must be doubled to 
obtain the actual turbine characteristics which have been 
plotted in Figures 3 through 8. 
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IV. DISCUSSION 



The agreement of both the predicted flow rate and the horse 
power obtained in Run 3 with turbine test data was encouraging 
It is to be noted however, that this agreement was obtained 
using a procedure which was conceptually incorrect. In Runs 
1 and 2, between-streamline continuity was checked and the 
streamlines were adjusted as necessary. In Run 3, between- 
streamline continuity was not checked, and as a result, the 
mass flow rate between streamlines was not precisely 251 of 
the total flowrate. It is noted however, that the deviations 
were less than 10.0% and while the radial positions of the 
streamlines varied by 10.%, the differences between predicted 
and measured output horsepower decreased from 24% to 4.5%. 
Since the enthalpy change on each streamline was computed 
using Euler's turbine equation, the total horsepower obtained 
by integration is sensitive to the streamline radial positions 
On the other hand, the calculation of the overall mass flow- 
rate is primarily a function of the blade throat openings and 
inlet conditions of the flow. Consequently, in relaxing the 
requirement for between-streamline- continuity , the output 
horsepower was changed significantly, while the overall 
flow rate was not. 
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Using this procedure, which preserves overall continuity, 
a performance map for the turbine was produced (Fig. 3-8) 
which agreed well with the off-design performance measurements 
made at lower speeds (Table III). It is noted however, that 
the inability of the program in its original form to predict 
the measured turbine performance is not explained, and both 
the program itself and the data input for the geometry should 
be closely re-examined. 

The difficulty in obtaining convergence to a solution at 
some operating points above the pressure ratio of 2.0 is 
likely to be the result of choking occurring on one or more 
of the streamlines. This was suspected but not fully explored. 

Finally, although the program was eventually made to run on 
the mini-computer, the time required to put the program into its 
final form was excessive since the original program was not 
written with segmentation in mind. When the segmented pro- 
gram was completed, only one operating point per run could be 
obtained. Thus, excessive time was spent compiling and loading 
the program. The execution time for the program averaged 2 
minutes at the lower pressure ratios and up to 30 minutes at 
the higher ones. This would be unacceptable if many points 
were to be examined. 
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V. CONCLUSIONS 



The program for calculating the performance of a single 
stage axial turbine reported by Macchi was revised, corrected 
and segmented and made to run on the Laboratory mini- computer . 
IVhen applied to the geometry of the air-drive turbine of the 
compressor test rig, selecting specific options for the repre- 
sentation of loss coefficients, the revised program failed to 
converge when design-point test conditions where input. Also, 
the computed horsepower was in error by as much as 24% when 
the program predictions were compared with specific test data 
obtained from the rig at off-design (lower speed) conditions. 
The revised program did however closely reproduce the results 
given by Macchi in his original report for a specific turbine 
geometry . 

When the requirement that the computed stream surfaces 
be such that they divided the flow exactly into equal 25% 
increments was removed, the program converged satisfactorily 
for design point conditions and gave agreement with test data 
to within 5% in flow rate and horesepower at off-design 
conditions . 

The complete performance map for the air drive turbine 
was obtained with the program following this revision. Based 
on the favorable comparison with data so far obtained, the 
map is likely to describe the performance to better than a 10% 
uncertainty. This is considered to be satisfactory for 



28 



sizing the turbocharger for the proposed compressor rig 
modification . 

The following recommendations are made concerning further 
application or development of the computer program: 

1. The failure of the program to converge before the 
final revision was made should be analysed closely, 
and the final revision removed if possible. 

2. The geometrical input for the air drive turbine (which 
was taken from drawings) should be reexamined and the 
physical dimensions of the blade rows themselves should 
be measured. 

3. Consideration should be given to putting the corrected 
original version of the program onto the IBM 370 com- 
puter so that, when successfully operating, a turbine 
map can be calculated with a single load. 
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TABLE I 



TURBINE GEOMETRY 

(see Figure 2; Dimensions in inches) 

STATOR: 





Hub Radius 


2.764 


Mean Radius 


3.196 


Tip Radius 


3.627 


Blade Chord 


1.003 


Blade Suction Side Radius of Curvature 


2..8065 



Maximum Blade Thickness .2252 



T.E. Projected Thicjcness .03 

T.E. Normal Thickness .0186 



ROTOR: 



Hub Radius 


2.693 


Mean Radius 


3.265 


Tip Radius 


3.837 


Blade Chord 


1.003 


Blade ''Suction Side Radius of Curvature 


2.8065 


Maximum Blade Thickness 


. 2252 


T.E. Projected Thickness 


.03 


T.E. Normal Thickness 


.0186 


Tip Clearance 


. 01 (estima 
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TABLE II 



MEASURED/DESIGN DATA USED TO VERIFY THE PROGRAM 



POINT 


RPM 


Tin 


(R)^OUT 


CR)^TO 


(PSI)P.R. 




H.P. 




1 


18310 


579.2 


550.8 


23.56 


1.602 


5.542 


110.1 


2 


15200 


557.4 


517.8 


20.43 


1.390 


4.698 


63.27 


3 


21300 


578.9 


506.8 


27.13 


1.846 


7.033 


172.0 


4* 


30500 


640.0 


— 


41.16 


2.8 


10.85 


485 



^Design Point 
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TABLE I I I 



COMPARISON OF PREDICTED TURBINE PERFORMANCE 
VS MEASURED PERFORMANCE 





FLOWRATE 




RUN I 


HORSEPOWER 




POINT 


PREDICT. 


MEAS. 


IDIFF. 


PREDICT. 


MEAS 


%DIFF. 




1 


5.88 


5.542 


6.09 


99.5 


110.1 


9.63 


2 


4.74 


4.698 


0.89 


52.5 


63.27 


16.17 


3 


7.04 


7.033 


0.009 


163.64 


172.0 


4.86 


4 


N.C. 


10.85 


— 


N.C. 


485 


— 


RUN 2 


1 


6.06 


5.542 


9.35 


90.92 


110.1 


17.4 


2 


4.90 


4.698 


4.29 


49.76 


63.27 


21.35 


3 


7.30 


7.033 


3.80 


130.76 


172.0 


23.97 


4 


N.C. 


10.85 


— 


N.C. 


485 


— 


RUN 3 


1 


5.82 


5.542 


5.01 


113.12 


110.1 


2.74 


2 


4.66 


4.698 


0.81 


61.96 


63.27 


2.09 


3 


7.04 


7.033 


0.10 


179.68 


172.0 


4.47 


4 


10.40 


10.85 


4.15 


444.18 


485 


8.42 



NC : Computer program would not converge to a solution after 

a large number of iterations. 
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TABLE IV 



VALUES OF ASSUMED TOTAL 


INLET TEMPEPA.TURE FOR EACH 


PRESSURE RATIO. GIVEN 


IN FIGS. 3, 5, 6, AND 7 


PRESSURE RATIO 


TOTAL INLET TEMPERATURE (°R) 




1.4 


545.5 


1.6 


562.6 


1.8 


577.3 


2.0 


591.0 


2.2 


603.6 


2.4 


615.3 


2.6 


626.1 


2.8 


636.6 



33 



ATMOSPHERE 




FIGURE 1: SCHEMATIC OF THE COMPRESSOR TEST RIG, WITH 

PROPOSED MODIFICATIONS 
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FIGURE 2: TURBINE ROTOR AND STATOR BLADE SHAPES 
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FIGURE 3: PREDICTED HORSEPOWER VS RPM AS A FUNCTION OF PRESSURE RATIO 

AT TEMPERATURES TO AVOID CONDENSATION 
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FIGURE 4: PREDICTED REFERRED HORSEPOWER VS REFERRED RPM AS A FUNCTION 

OF PRESSURE RATIO 
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FIGURE 6: PREDICTED HORSEPOWER VS PRESSURE RATIO AS A FUNCTION OF RPM 

AT TEMPERATURES TO AVOID CONDENSATION 
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FIGURE 7: PREDICTED FLOW RATE VS PRESSURE RATIO AS A FUNCTION OF 

RPM AT TEMPERATURES TO AVOID CONDENSATION. 
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FIGURE 8: PREDICTED REFERRED FLOWRATE VS REFERRED RPM AS A FUNCTION 

OF PRESSURE RATIO 



APPENDIX: A 



DESCRIPTION OF THE COMPUTER PROGRAM 
A-1. INTRODUCTION 

To enable the program to run on the laboratory computer, 
the program was divided into 4 parts; a main program and 3 
segments. A detailed discussion o£ program segmentation on 
the HP-1000 computer series is contained in Appendix C. In 
the description which follows, the program is treated as if 
it were one large program with many subroutines. 

The description follows the individual steps from 
start to finish in the analysis. A program flowchart is 
given in Figure A-1 and the FORTRAN symbols used in the 
program are listed in Tables A-I to A-IX. 

A- 2. DESCRIPTION 

A-2.1 Input Data 

There are 4 basic categories of input data; tur- 
bine geometry, operating conditions, special data and pro- 
gram control parameters. Since there was no card reader 
input device on the computer, all data were entered using 
either data or specification statements. Explanations of 
the turbine geometry, operating conditions, special data 
and program control parameters are found in Table A-I through 
A-V. The nomenclature for the blading is given in Figure 
A-2. 
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A- 2. 2 Initial Geometric Calculations 



The first calculation performed is to establish 
the 5 streamline locations at the stator inlet (station 0) . 
The streamlines are initially positioned such that there are 
equal areas (25% of the total flow area) between them. Next, 
blade heights of the stator and rotor are calculated using 
the hub and tip radii of each blade. Blade spacings for the 
stator and rotor are computed at 3 streamlines; hub, mean 
and tip. The blade spacing on the mean streamline for the 
stator is given by 

S = ^ Rm (A-1) 

s 

where S = Blade spacing 

Z = Number of stator blades 
s 

Rm = Mean stator radius 
A-2.3 Calculation of Gas Outlet Angles 

Subroutine VAVRA calculates gas outlet angles 
for both stator and rotor. The method is that of M.H. Vavra 
[Ref. 1] . The equation programmed in the subroutine is 






-1 




te 1 


^ 1 _ ( ^ ) \ 




s ^ 


S ' 


^0 / 



(A-2) 



where a = Gas outlet angle 

a = Throat opening 
S = Blade spacing 

tg = Projected trailing edge thickness 
This method is much simpler than that used by Macchi since 
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there is no variation in outlet angle with Mach ntiinber (for 
sub-sonic conditions). Therefore, once calculated, the stator 
and rotor exit angles remain unchanged. Subroutine VAVRA 
computes exit angles for the hub, mean and tip streamlines. 

The outlet angles at streamlines two and four are computed 
later in the subroutines STATR and R0T02. 

Before printing the input data, the program 
calculates the mean throat opening for the stator and for the 
rotor. The ten equally spaced radii and corresponding throat 
openings (part of the input geometry) are fitted with a fourth 
order Chebyschev polynomial. A throat opening corresponding 
to the mean radius is thus obtained. In the present appli- 
cation of the program to the drive turbine, the mean throat 
opening was obtained from the design drawing of the blading 
shown in Figure 2. It was assumed that the throat opening 
varied linearly with radial position and hence the throat 
openings at other radii could be calculated. The resulting 
throat openings are shown in the computer output under the 
heading of "Input Prints". The design values of the stator 
and rotor throat areas were obtained from the original 
design notes of M.H. Vavra. 

A-2.4 Calculation of the Flow Rate 

Subroutine CHAN is called to calculate the mass 
flow rate entering the stator. The equations used are as 
follows : 
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(A- 3) 



T- 

V ^ 

p - 

^ - 



T- 



TO 



I ^Mo*- 



-I 

F TO 



p/rt 

rr R Ti? ” R HozJ 



(A-4) 

CA-5) 
CA-6) 
(A- 7) 



Tf\ - p (A- 8) 

REF - J~l5Z (A-9) 

. «ro 1 g/c 

is the reference (dimensionless) flowrate 

and is used to check overall continuity later in the program. 

A- 2. 5 Solution of the Equation of Motion for the Stator 

Subroutine STAIR is called to solve the equation 

of motion for the stator outlet conditions. The equation of 

motion which is programmed is as follows : 



d CinY,’^) 

(A V.I 



" CO^ 




£r_\ 

L’- ' '' H u’- )' 



dS,* 

d^i 



2 Tf3ns < 



cjoi ' 

7^ 



s I I + 

y I 



Cl c\ M 

Y.^Vcj^ d^i 



C. Mcos^o^l ^ dSi^ 

— Sim 

Y| di-i (A- 10) 



45 




= « 
S 

! s 
s 
• 






where Ci “ J~ 



(a constant to convert H, the enthalpy 
from BTU 



IBM 



sec" 



Y, = 



X. -- 






V.^Cl') 

KCt) 



^ Axial velocity at a streamline 
Axial velocity at mean streamline 



Streamline radius 
Mean streamline radius 



fN m 


r 


1 - 






- ol 
dX 


Jbr\ 


Cl M cos^oi 1 




1 - 


Y. 

Cl Hco-s»o<(i-f ) 





C = Stator loss coefficient 

(which is initially assigned an 
estimated value) 

The derivation of this form of the equation of motion is 
given in Appendix B. However, at this stage of the analysis, 
the streamline curvature is assumed to be zero. Therefore, 
the equation of motion becomes: 



c^Xi 

dM . 



- - ^ I 



ciX. 






ciy. 



I _ C, H cos 

r. ^ 






ciXi 



C \ CoS I 



(A-11) 



The equation of motion is solved when the value 



of Y at each streamline satisfies the equation. The solu- 
1 

tion is to first put the equation in the form: 

<jUhY,“) 



dX 



= lU) 



(A-12) 



where I(X) consists of the right hand side of equation (A-11) 
Integrating equation (A-12) yields; 
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(A-13) 



r ^ 

- \ I W d y I ^ JU c 



where £.nc^ is the constant of integration when x = 1 and Y = 



1. With these boundary conditions Eq. (A-13) gives 

1 

Jbn -- - \ I Wdx. 



(A-14) 



yo 



using Eq. (A-14) in Eq. (A-13), 

X( 



j^y/= rCx)dX,-f xWclX 



which becomes 



r 

^Y, = xCxWX 



(A-15) 



(A-16) 



Taking the inverse natural log and the square root of both 

. X 

sides 



Y, = e 






X ex') d X I 



(A-17) 



Equation (A-17) is the form of the equation of motion solved 

in subroutine STATR. Solution of the equation gives five 

values of Y and thus the value of the axial velocity at 
1 

each of the five streamlines. Initially, the value of Y 

1 A 

is taken to be 1 and the value of is taken to be zero. 

1- 

In succeeding iterations, the calculated value of Y is used 

ds * 

to obtain a new value of and so on. 

1 

After calculating five values of Y^, the 'stator 
exit conditions are calculated at each streamline from the 
geometry of the velocity diagram. The convention for positive 
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and negative angles and velocities is defined in Figure A- 3, 
The required relations are the following: 



V<c, = Vo., C3) ■ Y, 


(A-18) 


\/xU| - \/cu i • T A Ny ^ 1 


(A-19) 


V, = \/cc, / cos o< « 


(A-20) 


Vg, = - Vo., [a R /.?L j 


(A-21) 



where L is the axial distance between stations and AR is the 



change in radial position of the streamline. radial 

component of velocity, is taken to be zero at this stage in 
the calculation. 

■fvTTv^ 



V, = 

T, - ^TO - 



V, 



■n.s - r, 



TO ~~ 



TC p 

TrO 






p. R. = P, / Ft 



P. - Pt( 



r... 



^■-1 



r, 



TO 



M, = V, /-I r<jc RT 



(A-22) 

(A-23) 

(A-24) 

(A-25) 

(A-26) 

(A-27) 



After the above quantities have been calculated at each 
streamline, subroutine STATR returns to the main program. 
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A-2.6 Calculation of the Stator Loss Coefficients 



The calculation of the stator loss coefficients 
at each streamline is accomplished by subroutine ALOSl. 

The method of solution to obtain these loss co- 
efficients is that formulated by Traupel [Ref. 7] . In 
Traupel's method, the value of the total loss coefficient is 
given by 



= C 









(A-28) 



total profile wall remaining 

The calculation of requires 9 subroutines. Figure A-4 

describes the connection between the subroutines and subroutine 



ALOSl. 



total profile 
to be 



The first step is to obtain the value of the 
loss coefficient, Cp • Cp is defined by Traupel 






(A-29) 



m 



m 



where C = initial value of the profile loss coefficient 
po ^ 

= mach number correction factor 
= trailing edge thickness correction factor 
= loss coefficient due to mixing losses and separa- 
tion losses 

= loss coefficient due to fan losses 

The total profile loss coefficient is calculated 
in the following manner. First, data for initial profile loss 
(Cp^) as a function of gas outlet angle (ai) for various values 
of gas inlet angle (ao) is read from an array (Fig. A-5). 
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This is done by subroutine TRAUl and functions XPO and YC . 



The values of are contained in two arrays XPOl (5, 8) and 



XP02 (6, 8). This is because the data shown plotted in Fig. 
A-5 has been divided into two sets. One set is for values 
of Ui between 40° and 80°. The other is for values of 
between 80° and 170°. The FORTRAN symbols for the two ranges 
of values of oti are ALFOl(I) and ALF02(I) respectively. The 
FORTRAN symbol for the gas inlet angle is ALFl (J) once the 
data pints selected from the plot are entered, fifth and sixth 
degree polynomials respectively are fitted through the data 
points. The value of can then be determined for given 
values of ai and Uq. 

The mach number correction, X is obtained from 

m 

Fig. A-5. Subroutine CSIM calculates the value of X^ using 
straight line approximations of the plot. 



tained from the data in Fig. A-6 using the linear interpolation. 

The abscissa of the curves for X and I is either f or 1-f 

s m 

where f is defined as 



where 6 = normal trailing edge thickness, 
t = blade spacing. 
ai = gas outlet angle. 

The loss coefficient due to wall friction, C , 

w 

is calculated using 




Subroutine CID calculates the remaining terms in 

n for C • These are X , C , Cr* They are ob- 
^p s m f ^ 





(A-30) 
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CA-31) 





~t SlKJoC 



where t - blade maximum thickness 
= blade height 

This equation is programmed in subroutine CSIW. 

The value of is obtained using subroutine CSIR. 
Cr is defined by Traupel to be an all-inclusive loss coeffi- 
cient which accounts for any remaining losses not previously 
defined. It is written as 






RO 



(A-32) 



is an initial value of which depends on the value of 
<p, where ({> is given by 



(|) = 



Vi SllOo< 



^J (A-33) 

in which Vi = true velocity of gas 
V = blade speed 

A plot of vs Ip is shown in Fig. A- 7. The correction 
is a function of s/Z where 
s = chord length 
Z = blade height 

and is obtained using the data in the lower half of Fig. A- 7. 

The total stator loss coefficient is computed 
for 3 streamlines; those at the hub, mean and tip. 

The loss coefficients at streamlines 2 and 4 are 
obtained by linear interpolation. 
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A refinement to the stator loss coefficient may- 
be applied depending on the input value of one program control 
parameter. The following 3 variations of 5^ are available: 
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(A-35) 



(A-36) 



where C 

0 



loss coefficient calculated using the method of 
Traupel Y-l 





(A-37) 



The values of the program control paramenter 
required to select between options are given in Table A-V. 
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* 



I- 







I 





Before returning to the main program, subroutine 
ALOSl calculates a value of 5* which is a blockage factor to 
be used in the equation of continuity. There are three ways 



to define ?*; they are as follows: 




(A-38) 



(A-39) 




(A-40) 



A-2.7 Solution of the Continuity Equation After Re- 
turning to the Main Program 

The overall continuity at the stator exit is 
checked. Subroutine FLOWR performs this task. The flow chart 
for FLOWR is given in Fig. A-8. In FLOWR the mass flowrate 
required by continuity is checked against the calculated mass 
flowrate. If the calculated flow rate does not agree with 
that required by continuity, adjustments are made to the axial 
velocity and/or the inlet Mach number, as will be explained. 



The mass flow required by continuity is 



REGlI) 



REP 




(A-41) 



where = reference mass flovrrate as computed in subroutine 



CHAN 
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= # of stator blades 
= mean stator throat opening 
= mean stator radius 

The mass flow rate at each streamline com- 
puted in this subroutine is 




CA-42) 



where Z is an area reduction coefficient defined by 




H***- I 




(A-43) 



Z gives the percentage of flow area between the blades over- 
which it is permissable to assume a uniform velocity. The 
boundary layer on both sides of the flow limits the available 
flow area and the backage facotr, Z. accounts for this. 
Equation A-43, Z is seen to be a function of the energy para- 
meter H*** and C*. 5* is the value of the loss coefficient 

returned from subroutine ALOSl. The energy parameter is 
defined as 
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u*** - = Energy thickness 

~ Displacement thickness 



(A-44) 



H*** can be written as 



= L. 
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where ; 



and 



m = .15 
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for unchoked flow 



for choked flow 



The derivation of Z and H*** os given in Appendix B. 

The expression for $ , the flow function, 

for unchoked flow is 
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(A-46) 



and for choked flow is 
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(A-47) 
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After calculating for each streamline, the flow 
rate is integrated from hub to tip and the resulting value is 
compared with If the two values of flow rate agree to 

within a specified tolerance (see Table A-IV) continuity is 
considered to be satisfied. Then, after calculating the total 
percentage of mass flow between adjacent streamlines, sub- 
routine FLOWR returns to the main program. 

If the flow rates are not within tolerance the 
program checks to see if the actual mass flow is to high. If 
it is to high, the value of the axial velocity is lowered 
proportionally to the difference between the actual and 
required flow rates . 

If the actual flow rate is too low, the procedure 
is more complicated. First, the flow is checked to determine 
whether choking has occurred. Streamlines one and five are 
checked. If the flow is in fact choked at those streamlines, 
the inlet Mach number is lowered and the program loops back 
to recompute the reference mass flow rate and repeat the 
complete procedure. 

If the flow is not choked, the axial velocity is 
raised proportionally to the difference between actual and 
required flow rates and subroutine FLOWR returns to the main 
program. 

A-2.8 Calculation of the Rotor Inlet Conditions 

Continuity having been satisfied through the 
stator, the rotor relative inlet conditions are calculated. 

In subroutine ROTOl, the following expressions are used: 



56 



u = 



gj R 
I ^ 



= 



UJ 

12. 



R 



STfVroR. 



R 



RoTOft. 



WUu, - Vju.^ -u 



/3 



I 




Wm., 

Vo.. 










( r. * w, ) 

\ cp / 






r-i 



(A-48) 



CA-49) 



(A-50) 



CA-51) 



(A-52) 



(A-53) 



(A-54) 



(A-55) 



57 



(A-56) 



He “ 



Inhere and are equivalent temperature, pressure 

and enthalpy respectively. 

A- 2. 9 Calculation of the Rotor Exit Conditions 

Calculation o£ the rotor exit properties follows 
the same procedure as was used to compute the stator exit 
properties. The process is outlined here with notable 
differences explained. Subroutine R0T02 calculates the 
rotor exit properties. A flowchart of R0T02 is given in 
Fig. A-9. 



of motion for each streamline. The equation of motion in 
terms of relative quantities is 



The first step in R0T02 is to solve the equation 






Cl He ol 



(A-57) 
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At this point in the calculation, streamline curvature is 
neglected. Hence, Eq. (A-57) reduces to 



The derivation of Eq. (A-57) is contained in Appendix B. 

Equation (A- 5 5) is similar in form to Eq. (A-10) . Hence, the 

method of solution is identical to that employed by subroutine 

STATR. However, after solving the equation, the value of 

at each streamline is examined to determine whether or not it 

falls into the range .2<Y <2.0. Values of Y greater than 

2 ^ 

2.0 are set equal to 2.0 while those less than .2 are set 

equal to .2. Successive values of Y^ at each streamline are 

compared, and when the values of successive iterations are 

within a specified tolerance (see Table A-IV) , the iteration 

ends . The values of Y are used to calculate the rotor exit 

2 

conditions using the following equations: 
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(A-60) 
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Subroutine R0T02 , then returns to the main pro- 



gram. 



After calculating the rotor outlet conditions, 
the rotor loss coefficients are computed. Subroutine AL0S2 
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calculates the rotor loss coefficients following the process 
used in subroutine ALOSl for the stator losses. The principle 
exception is that a tip clearance loss is also calculated and 
added to the total loss coefficient. The tip clearance loss 
coefficient is obtained from subroutine ALEAK which uses a 
straight line approximation to the curve shown in Fig. A-10. 
Subroutine AL0S2 also computes values of and one of the 
three refinements to 

Subroutine FLOWR is called to check continuity 
at the rotor exit. If continuity is satisfied, the program 
continues. If not, the same procedure is followed as pre- 
viously explained for the stator outlet (Fig. A-1) . 

A-2.10 Accounting for Streamline Curvature 

All calculations to this point have neglected 
streamline curvature and assumed that the streamlines remain 
fixed through the stator and rotor (Fig. A-11). The radial 
shift in a streamline between stator inlet and rotor outlet 
can be written as 

R •trrfvToR. ^ RoToR (A-67) 

IKJUET (XTTLET 

This is the net radial shift in a streamline between stations 
'O' (stator inlet) and '2' (rotor outlet). It is shown in 
Section 16.4 of Ref. [5] that the radial shift in a streamline 
between the stator and the rotor (station 1) can be written 
as 
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strtoR 

OOTue.T 
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VTRTDfi 

JKJLST 
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ovTlBT 



(A-68) 



The angle between the meridional velocity V and the axial 
° ^ m 

velocity is X. Tlte radial velocity can be expressed as 



= Vo. X 

and from Fig. 16(1) of Ref . [5] , it follows that 



(A-69) 



- A R 



Tam X = a^i_ 



(A-70) 



AR 



Using Eq. (A-68) in Eq. (A-67) , 

“ Vo, 

where = Average streamline slope 

Also, from using Eq. (A-68) 



a L 



(A-71) 



Rearranging ; 



Cos X = 



21 . 





(A-72) 



COS X 




(A-73) 
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The remaining term used in the calculation of streamline cur 
vature (Section 16-4 of Ref. [5]) is 



K 



6R 

TJ 



where K is the so called curvature factor. It usually has 
a value between 4 and 6 and in the program its value is taken 
to be 5. Having calculated cosA, AR and 6R, the program 
repeats the solution process. However, the only quantity 
which is unchanged is the reference mass flow rate In 

subroutine STATR the equation of motion is solved, this time 
accounting for streamline curvature. The same is true in 
subroutine R0T02. 

The flow path of the program is identical to the 
section which did not account for streamline curvature. Next 
the program computes an average pressure ratio at the rotor 
outlet using the expression 



^TO 



-] -- 
fro /STReAMUKiE 4 





Si- 

2 



(A-74) 



If this pressure ratio is within a specified tolerance to 
the actual pressure ratio (which is input data) the program 
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proceeds to the final stage of the calculations. If not, the 
inlet mach number is adjusted by an amount which depends on 
the difference between the calculated and specified pressure 
ratios. If the calculated pressure ratio is too high, the 
Mach number is lowered using 



If the computed pressure ratio is too low, the Mach number 
is raised using 



In both cases, the program loops back to subroutine CHAN and 
proceeds to compute a new reference mass flow rate based on 
the new value of the inlet Mach number. The entire process 
is then repeated until the pressure ratios agree within the 
specified tolerance. 

A-2.11 Final Calculations 




Pressure Ratio Difference 
18 



(A-75) 




Pressure Ratio Difference 

n 



(A-76) 



Stator and rotor outlet conditions not previously 
calculated are computed as follows: 



o< a. - Tan/ 







(A-77) 




Vo.1. 



(A-78) 



COS <^ 2 . 
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(A-80) 



A =• 




ROTOR EXIT 
RELATIVE MACH # 



W<3. 



CA-81) 



(A-82) 



(A-83) 



(A-84) 
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d. - 


/ To ' Itt, 


(A-88) 


( r-s 


\ 

o 


Stator Blade 
Efficiency 


1 1 


(A-89) 


Rotor Blade 
Efficiency 
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(A-90) 
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(A-91) 


Head Coefficient 


;2 9c T C'Pro -Tai<>^ 


(A-92) 




Blade-Jet Ratio 


-1 

= [Head Coefficient] 


(A-93) 



Stator Exit 
Relative Mach # 



W> 

- YRg<iTi 



(A-94) 
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The turbine horsepower is obtained by integration. 
The Ah term at each streamline is weighted by the percentage 
of mass flow at that streamline. The product is then inte- 
grated from hub to tip and result, Ah, is used in the tur- 
bine horsepower equation 

« 

A k ’X' 

5TO (A-95) 



H.p. = 



The moment is calculated using 



M = 



( H.p. X fro) 



oo 



(A-96) 



Referred horsepower, moment, mass flow and RPM are calculated 
using 



H.P. 



(A-97) 



I 

^REF - 

RPM 



M 

m 0 

5 

RPM 

0 



(A-98) 



(A-99) 



(A-lOO) 
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where 



© 



T 



TO 






Pro 

IH.7 

The values of the total-static efficiency, total- 
total efficiency, total-static pressure ratio, total-total 
pressure ratio, head coefficient, blade/jet ratio, r* and 
inlet mach number are then averaged. 

With all calculations completed, the results are 
printed under the heading "STATOR SOLUTION”, "ROTOR SOLUTION", 
and "OVERALL TURBINE CHARACTERISTICS". 
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TABLE A- I 





TURBINE GEOMETRIC INPUT DATA (STATOR) 




(see Figure . 


^^-2; Dimensions in inches) 


FORTRAN 


SYMBOL 


DESCRIPTION 




ZS 




Number of blades 


RSCD 




Hub radius at stator outlet 


RSC3) 




Mean radius at stator outlet 


RS(5) 




Tip radius at stator outlet 


C 




Blade chord (mean) 


Cl 




Blade chord (hub) 


CO 




Blade chord (tip) 


E 




Blade curvature (mean) 


El 




Blade curvature (hub) 


EO 




Blade curvature (tip) 


T 




Maximum blade thickness (mean) 


TI 




Maximum blade thickness (hub) 


TO 




Maximum blade thickness (tip) 


TE 




Projected T.E. thickness (mean) 


TEI 




Projected T.E. thickness (hub) 


TEO 




Projected T.E. thickness (tip) 


TN 




Normal T.E. thickness (mean) 


TNI 




Normal T.E. thickness (hub) 


TNO 




Normal T.E. thickness (tip) 


Al(l-lO) 




Ten values of throat diameter at 
10 equally spaced radii 
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FORTRAN SYMBOL DESCRIPTION 





AL 


Blade camber line length (mean) 


ALI 


Blade camber line length (hub) 


ALO 


Blade camber line length (tip) 


RC(1) 


Hub radius at stator inlet 


RC(3) 


Mean radius at stator inlet 


RC(5) 


Tip radius at stator inlet 
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TABLE A- I I 





TURBINE GEOMETRIC INPUT DATA (ROTOR) 




[see 


Figure A-2; Dimensions in inches) 


FORTRAN 


SYMBOL 


DESCRIPTION 




ZR 




Number of blades 


RRCD 




Hub radius 


RRC3) 




Mean radius 


RRC5) 




Tip radius 


CR 




Blade chord (mean) 


CIR 




Blade chord (hub) 


COR 




Blade chord (tip) 


ER 




Blade curvature (mean) 


EIR 




Blade curvature (hub) 


EOR 




Blade curvature (tip) 


TR 




Maximum blade thickness (mean) 


TIR 




Maximum blade thickness (hub) 


TOR 




Maximum blade thickness (tip) 


TER 




Projected T.E. thickness (mean) 


TEIR 




Projected T.E. thickness (hub) 


TEOR 




Projected T.E. thickness (tip) 


TNR 




Normal T.E. thickness (mean) 


TNIR 




Normal T.E. thickness (hub) 


TNOR 




Normal T.E. thickness (tip) 


TIPC 




Tip clearance 
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FORTRAN SYMBOL 


DESCRIPTION 






A2(l-10) 


10 values of throat diameter at 
10 equally spaced radii 


ALR 


Blade 


camber line length 


(mean) 


ALIR 


Blade 


camber line length 


(hub) 


ALOR 


Blade 


camber line length 


(tip) 


CV 


Axial 


distance between stations 


CK 


Curvature Factor 
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TABLE A- I II 



FORTRAN 


TURBINE OPERATING CONDITIONS (INPUT DATA) 
SYMBOL DESCRIPTION 



AMC 


Assumed inlet Mach number 


AMS 


Assumed stator exeit Mach number 
(absolute) 


AMR 


Assumed stator exit Mach number 
(relative) 


PTO 


p 

Total inlet pressure ( TO) 


TTO 


T 

Total inlet temperature ( TO) 


PR 


Total-static pressure ratio 


RPM 


Operating speed (RPM) 


VA1(3) 


Assumed axial velocity in stator 


VA2C3) 


Assumed axial velocity in rotor 
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TABLE A- IV 



SPECIAL INPUT DATA 



FORTRAN SYMBOL DESCRIPTION 





TOL 


1 


Toler aice for convergence of 
equation of continuity 


TOL 


2 


Tolerance for between-S.L. 
continuity (not used) 


TOL 




Tolerance in pressure ratio 




convergence 


TOL 


A 


Tolerance in equation of motion 


H* 


convergence 
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TABLE A-V 





PROGRAM CONTROL 


PARAMETERS 


FORTRAN SYMBOL 


POSSIBLE VALUE 


EFFECT/MEANING 




IND 


1 


Prints results in sub- 
routines STATR, FLOIVR., 
ROTO 2 




1 


No printing in the above 


ICL 


1 


Rotor is shrouded 




1 


Rotor not shrouded 


ICOZ 


1 


^ = Co 




6 


^ C (Y Pressure Ratio) 




8 




ICON 


1 






2 


^ “ ^PROFILE 




3 


^ ^ ^TOTAL 
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TABLE A-VI 





FORTRAN 


SYMBOLS IN THE MAIN PROGRAM 


FORTRAN 


SYMBOLS 


DESCRIPTION 




BESP 




3* = [1 + ^ 


01 




Stator throat opening (hub) 


00 




Stator throat opening (tip) 


OIR 




Rotor throat opening (hub) 


OOR 




Rotor throat opening (tip) 


0 




Stator throat opening (mean) 


OR 




Rotor throat opening (mean) 


ANG2I 




Stator gas outlet angle (hub) 


ANG20 




Stator gas outlet angle (tip) 


BETAI 




Rotor gas outlet angle (hub) 


BETAZ 




Rotor gas outlet angle (tip) 


G 




Gray, constant, 32.174 FT • 2 

\ml ^ • S 6 C 


CJ 




778.16 FT.LBF/BTU 


EXPl 




r—H 

1 


EXP 2 




T-Vy 


ERRE 




Gas constant, 53.3459 

^ LBM.oR 


EMME 




Molecular mass, 28.970 LBM/LB MOLE 


GAM 




Y, Ratio o£ specific heats 


ETAT 




Total-total efficiency 


ETAI 




Total-static efficiency 


ETAS 




Stator blade efficiency 
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FORTRAN SYMBOL 


DESCRIPTION 




ETAR 


Rotor blade efficiency 


RSTAR 


Theoretical degree of reaction 


ALOS 


Head coefficient 


BLAJE 


Blade/jet ratio 


DRl 


Radial shift of steamlines 


AMWl 


Stator exit relative Mach Number 


AMSl 


Stator exit absolute Mach Number 


AMV2 


Rotor exit absolute Mach Number 


AMR2 


Rotor exit relative Mach Number 


DELH 


A 


HP 


Horsepower 


AMOM 


Moment 


THETA 


9 


DELTA 


6 


HPl 


Referred H.P. 


AMOMl 


Referred moment 


RPMl 


Referred RPM 


WLBMl 


Referred mass flow rate 


ETAS 


Average total-static efficiency 


BETA6 


Average total-total pressure ratio 


ETA6 


Average total-total efficiency 


AKIS5 


Average head coefficient 


RSTAR 5 


Average theoretical degree of reaction 
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I 





tc 






TABLE A- VI I 



FORTRAN 


FORTRAN 


SYMBOLS IN SUBROUTINE CH/\N 


SYMBOLS 


DESCRIPTION 




TTO 




T^q, total temp, at station 0 


AMC 




Inlet Mach number 


PTO 




P^Q, total pressure at station 0 


RC (I) 




Streamline radii at station 0 


WLBM 




• 

M, required mass flow, pAV 


TC 




Static temperature 


VC 




Velocity 


PC 




Static pressure 


RHO 




p, density of air 


WCHAN 




Mref> reference mass flow 


WPERO 




% of M at each streamline 
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TABLE A-VIII 



FORTRAN 


FORTRAN 


SYMBOLS IN SUBROUTINE STATR 


SYMBOL 


DE SCRIPTION 




ALFAl 




Stator gas outlet angle 


Y 




Ratio o± streamline radius/ mean 


A 




radius 


AMS 




Mach Number at station 1 


T 




Static temperature 


P 




Static pressure 


VI 




Absolute velocity 


VAl 




Axial velocity 


Y 




Ratio of axial velocity to mean 






axial velocity 


s 




Entropy 


DSDX 




Entropy gradient between streamlines 


VUl 




Tangential velocity 


PRAT 




(Total-static pressure ratio) ^ 


TIIS 




^IIS 


DALF 




da 

3x 


RSF 




Mean stator radius 


DELR 




D p 

Stator in - 'rotor out 


ZETAPS 






ZETAS 




§ 

s 


VRl 




Radial velocity 





79 



TABLE A-IX 



FORTRAN SYMBOLS IN SUBROUTINE TRAU2 



FORTRAN SYMBOL 


DESCRIPTION 






GSIP 


Xp , correction to 


R 


initial profile loss coefficient 
po ^ 


ANGl 


Gas outlet angle 


ANG2 


Gas inlet angle 


R1 


X , Mach No, correction 
m 


R3 


§p, remaining loss coefficient 


R2 


g loss coefficient due to wail 

^W’ friction 


RPRO 


§p, total profile loss coefficient 


CL 


Rotor tip clearance 


YCL 


Tip clearance loss coefficient 


RTOT 


Total loss coefficient 


T 


Blade spacing 


DEZ 


Normal trailing edge thickness 


HM 


Blade height 


CSID 


„ trailing edge thickness 

■^6’ correction factor 


PSID 


§^, loss coefficient due to fan 




losses 


PSIF 


§ , loss coefficient due to mixing 
^ and separation- . 


UM 


Tip speed 
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TABLE A-X 



FORTRAN 


SY^ffiOL SIN SUBROUTI NE FLOWR 


FORTRAN SYMBOL 


DESCRIPTION 




PRATCR 


Critical pressure ratio 


PHICR 


$CRiT» critical flow function 


HSTAR 


H***, energy parameter 


XI 


Z, area reduction coefficient 


PHI 


$, flow function (unchoked flow) 


ARAT 


Streamline throat DIA/mean throat DIA 


SUM 1,2, 3, 4 


Successive values of the 
flow integral 


AS 


Mean stator throat diameter 


AR 


Mean rotor throat diameter 


WREQ 


• 

M required to satisfy continuity 


WSUM 


M calculated 


WPER 


% of M at each streamline 
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TABLE A-n 





FORTRAN 


SYMBOLS IN SUBROUTINE ROTOl 


FORTRAN 


SYMBOL 


DESCRIPTION 




OMEG. 




w, wheel speed (PJVD/sec) 


U 




wR . . 

stator mean 


U2 




R * 

rotor mean 

^ R 

* stator mean 


WUl 




Wpi see figure A-3 


BETAl 




6 1 , see figure A- 3 


W1 




Wi, see figure A-3 


TTE 




Equivalent temperature 


PTE 




Equivalent pressure 


HE 




Equivalent enthalpy 


ZETAR 




rotor loss coefficient 


ZETAPR 




§ , rotor profile loss coefficient 
P 


DHEDX 




Enthalpy gradient between streamlines 


DSDX 




Entropy gradient between streamlines 
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TABLE A- XI I 



FORTRAN 


FORTRAN 


SYMBOLS IN SUBROUTINE R0T02 


SYMBOL 


DESCRIPTION 




BETA2 




B 2 > see figure A-3 


DBETDX 




^ between adjacent streamlines 


VA2 




V , axial velocity 
a. 2 


W2 




W 2 , see figure A- 3 


CL 




Axial distance between stations 


WR2 




Radial component of velocity 


WU2 




Wu 2 , see figure A-3 


VU2 




Vu 2 , see figure A-3 


AMR 




Relative Mach No. at rotor exit 


T2 




T 2 


T2S 




T 






2S 


P2 




P 2 


PPJ^T2 




[Total-Static pressure ratio] ^ 
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FIGURE A-1: PROGRAM FLOWCHART 
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0 « THROAT DIAMETER 
E » CURVATURE RADIUS 
S » SPACING 
T » MAXIMUM THICKNESS 

TE « TRAILING EDGE THICKNESS PROJECTED IN 
PERIPHERAL DIRECTION 

TN « TRAILING EDGE THICKNESS, NORMALTO FLOW 
DIRECTION 



FIGURE A-2: BLADE NOMENCLATURE 
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FIGURE A- 3: VELOCITY DIAGRAM NOMENCLATURE 
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MAIN PROGRAM 




FIGURE A-4: INTERCONNECTION OF THE SUBROUTINES IN THE 

TRAUPEL METHOD 
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AtA- g.42 G nm d w crt Cpo ^ Pnrf2verluatei for BcschlcTmigimgsgitter raid Machzahlkorrektnr Xjj 



FIGURE A- 5: INITIAL PROFILE LOSS COEFFICIENT AND MACH NUMBER 

CORRECTION FROM TRAUPEL 
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Abb. 8.4.4. Korrekturfaktor X9 MLjchverlost Cm 
ixLfolge endlicber Austxittskantendicke Oder 
Ablosung 




FIGURE A-6: T.E .THICKNESS CORRECTION FACTOR, MIXING LOSS 

COEFFICIENT AND FAN LOSS COEFFICIENT FROM 
TRAUPEL 
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FIGURE A 



(jntere Werte nurbei Profile/), 
ale wenig empfinwich 
oaf Zustrdmwinkelanderung 
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Abb. 8.4.7 Rajidverlust in Tnrbinenschaufelungen 



-7: ’’REMAINING" LOSS COEFFICIENT FROM TRAUPEL 
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FIGURE A-8: SUBROUTINE FLOWR FLOWCHART 
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FIGURE A-9: 



SUBROUTINE ROTO 2 FLOWCHART 
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FIGURE A-10: TIP LEAKAGE LOSS COEFFICIENT PLOT FROM TRAUPEL 



I 
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FIGURE A-11; STREAMLINE COORDINATES 



APPENDIX: B 



DERIVATION OF EQUATIONS USED IN THE PROGRAM 
B-1. EQUATION OF MOTION FOR RELATIVE FLOW: 

The equation of motion for relative flow Ref. [5] is 

VHr * W X (v X W ->-o2co) h TVS 

(B-l) 

Using cylindrical coordinates, the terms of EQN (B-l) may be 
expressed as follows : 



- L& 



(^0 ^ Ht?. , / 
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^ r 

► 





+ 



C»r 






(B-3) 
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(B-5) 
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TVS = T 



C(5 . 






r 9o 



■ ^ • ?s 

+ Cj + Cf 



(B-6) 



Combining equations (B-l) through (B-6) the terms in (B-2) 
can be written as: 



_L ^"^ 0 ^ _ ;^(hwJ 

r a© ' r a© ?2- 

Wr 

■ - 1 — II. ■ • 

r 

(B-7) 



d (v-WxQ 

^)r 



^W, 






ir 9© 



5Mb. 

92r 




32- 
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9(r WxJ 
9^ 
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7- ^ 
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9 Hr „ ^(^WjLc) 
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(B-8) 
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(B-9) 
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Since the flow has been assumed to be axisymmetric, all deri- 
vatives with respect to 0 are zero. Thus, Equations (B-7) , 
(B-8) and (B-9) reduce to, respectively: 



o = 



Wo. ^ (jr Wxu) 

T 



Wf > (rWxu) 
r h'r 



a LU Wv- 
(B-IO) 



9 Hr _ Yv/ _ v\/ 

^^2: ■ ^2: ^ r 

9 Wav] ^ 

‘ ' 9i- 



(B-ll) 



^ H R ^ V\/u^ ^ Wr ^ ^\l\Jou , 

r ^ 5t" 

2ouW^i-T 

(B-i 



Equation (B- 


10) may be 


written as 








_ __ Wr 


3. 


3. 


Wy 


9 ^ 


Wo, 


9r 




Wco (B-13) 


Substituting 


into equation (B-ll) , 




- 


9Hr 




Wr JUWau) 


■ f VVr 


9\A/r 




r 


WCK, 9>" 




9 


Wf 


- ;au ■ 


Wu^Wr -p 


9S 




9r 




Wc^ 


32- 


(B-14) 
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Multiplying equation (B-9) by Wr and Equation (B-14) by Wa 
results in 



^ ? Hr _ WxcWr ^(.fWxJ) \A/. Wv j. 

*■ ■ I- 



Wa.Wr + a UJ 4- Wf r ^ 



and 



Wo 






3Hg. — 






(B-15) 
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(B-16) 



Adding these two equations yields 



9M( 



Wr + Wo. - r 






^Ir 



5 ^ 



\A/r -^+^0, 

(9^* 3 i 



CB-17) 

Since the flow has been assumed to be adiabatic, the total 

relative enthalpy, H , is constant along a streamline. Thus, 

K 



vHr 

and re-arranging. 



<>' CB-18) 



Wr d Hr 



32 
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From equation (B-19), eq. (B-17) can be written as 
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32 
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(B-19) 



(B-20) 
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Substituting Eq.s (B-19) and (B-20) into equation (B-15) gives 



Wr 


Wxc V\/r 
h Wcx, >r 


). w 


- 




V\/A-u^A/y- 


W»- -r 








Wo. 

-w 


Wo, 




Multiplying Equation 


(B-21) by -j-2. 

R 


yields 




PHr _ 


Woe 


H^-\A4u.') 


- W^- 


5 V\/r 


av- 


V- 










(B-21) 



-i- 



Wou ^ ;luo Wxo i- r 



21 



(B-22) 



This expression is identical to equation (B-21) and is the 
equation which must be solved. It must be put into a form 
which can be solved by the computer. Re-writing equation 
(B-22) given that 






cc 



Wcx. _ _L ^ C.wJ’) 

^ ^ |r 



yields 



3Wr , 2W^ jCir'WjJl , 

-57 2W^-^ + — p ^ 

H U> + ST 

(B-23) 



The relative enthalpy can be written 
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(B-24) 
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The equivalent enthalpy, defined in ref. [1] is 

u l . ^ ^ 

Hence, the relative enthalpy can be written as 



CB-25) 



Hr - He - 



2. 

Ua 






(B-26) 

Also, the turbine outlet static temperature can be written 



as 



r, - 



H, 



W; 



Cp cS.C'p 



(B-27) 



Substituting Eq. (B-26) and Eq. (B-27) into Eq. (B-23) and 
applying Eq. (B-21) to the rotor exit, gives 



M cvWiLi-i -i. 



Wa‘ 



5Cp 





n 






^ J 


+ :i 



Ml 






9S; 



3tr, 



= O 



(B-28) 



Given the relationships: 



r«M = A = 






101 



and 



1 4- - 



_J 

Cos^'X 



Equation (B-28) can be written as 



<? Wo. — Wg,x ^ 

^ ^ ^ Cp cos^ X 2. ^ 






-f- 4 cu vvCic , - 1 

iTx >r-j_ * ;ir2_ 



+ 



ACL!i) 






iHs ~ V\/^ 



and substituting 



ii^=3Co\ i 






(B-29) 



^ na 



gives 



^ CWa,!! 
5r, 



into equation (B-29) 






- Wc^, . 

^ Cp cos*Xt- 



.p HcuWA>ua.-2 









-P 



a • 

a CU T;l ^- 



'? L 












(B-30) 



r 

Multiplying Eq. (B-30) by (-^-^) results in the dimension- 

% 



less form of Equation (B-29): 
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lfr\ 



w, ^ 






Am 






- a 



7 ( WrJ 

Wcu*, V\/ci,^ 






V^o3a 



Pz, WALa.W<c^ 
^ ^ WgA.A Wo^a, rznv ^ ^ L^^mWa,^^Wcc^ 



Vy/oJ^ CpCo^l^^i Wia^^^Wo^ V\ ^ 

•fA ' 

^n turamW-^tu^Wq.^ _ ^ r^ry^ ;^Hg 

\A/.^ Wft., an. 



2.Wj n. Ifarn , 

" w4^ " 



rjtm 



a. Me 

WcJ 



Am 



W>/g>a W^oua 

Wq!, 



ISx 



CB-31) 



Introducing the non-dimensional quantities 



Y = 



W< 



A/ 



Wco 



m 



CB-32) 



y - 



n 



yv\ 



CB-33) 







(B-34) 
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Equation (^-31) is written as 



9CY^) , ^ 

• “ ot " — : • rr\ 









Wc 






co5.^ 'x ay 



+ 



Tft^/3 a CxY T-ftN<g ) L, U^YTptog 

X ax ^ ^ 



w, 



aHe 4 . 



ou 






ay 






as 



ax 



^4, 



= O 



avte 



(B-35) 



The fourth term of Eq. (B-35) is 



^ ( XYTftM/g) _ 

Y- d'Y ^ 



X Y 



31^ + y TflN ,6 4^- ^ Y rA^J ^ 4^ 



ax 






;iY*rpv'j§ 1- ^4. 4 - 
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also , 



j Tp>n 3 _ 



ay 



■"'zyT.»'i3 



Cos (S 



- Tan T 






jCy') 

ay 



Therefore, equation CB-35) can be written 



^y ^ ‘ ^ Wo. "" n COS' 



T TftN 3 a^ , ^ 

'^' cos*(3 av y 






cos^'v ay 

4U,KYrftM& 



+ 



vy<v 



nr\ 



w/ 9-y 



91-Ie. ^ Ur^VJ-i. 

•h 2. — + 









;?He _ 






3s 



ay 



(B-36) 



Multiplying Eq. (B-36) by ( ' ^ observing that 

( I + TP»M* (S = ) j ^ 



J ^CY") ^ . 2. 
a-y ^ 



5^ rf»% ^Vv V- 



Wk 



I 






<^2r cos^x a-y 



2Tp,io^-^-h ^ 






y 



^ ^ COS 8 

• \ K / \^y 1ft./ ^ 









a.cos^^ ane 



vv, 






s^y 



+■ 



2 M 6 CoS 3 

Y= 



-Si 



21_ - O 

av ^ 



(B- 37) 
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To account for streamline curvature the following terms are 
introduced: 




L 



2 . 








5. 



(B-38) 

where X, the angle between the axial and radial components of 
velocity at a point, is approximated as the average value be- 
tween two stations. 

Also , 

K S R _ _ » 

L* ?• 2" 

(B-39) 

where 6r is the streamline shift throught the rotor defined 

^ ^ROTOR - ^ROTOR 

OUTLET INLET (B-40) 

Substituting Eqs . (B-38) and (B-39) into (B-37) yields 









” X Tan (3 



1 

^ « 
Lf 

1 1_ 


' L^+ CARf ) 


'' L-‘ / \ 


L* / dy 



d 3 
dX 



_2 

X 



^ ^ t KJ CoS ^ 






;2Ufr^Ua cos^/3 P cos d He 

Wo,’Y^ 'it 



<3 He cos*& 

W/ Y “ 

P\ 



Sin’^ 



A" 

d)C 



(B-41) 



106 



To obtain a dimensionless equation, the term 




is introduced into Eq. (B-41) giving 








^ CSO Cm \ 




1 ci %* 


)' 


j 


' dy 



-,-P a dp 



3. _ Hld.K^i»o|Sco5g 



3.Vr^V3.Cl/i ^ 

^ m ^ 



Ct cos^^ d N& 

dY 



Ct Me CoS ^ 




ds! 

dX 



(B-42) 

Equation (B-42) is the form o£ equation of motion used in the 
computer program. 



B-2. EQUATION OF MOTION FOR ABSOLUTE FLOW 

The equation of motion for absolute flow 



= y (vxv)-i-TVS 



(B-43) 
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Differs from the equation of motion for relative flow 



y -t-54Aj) + TVS 

CB-44) 

only by the term i^x2w which is the Coroiolis acceleration. 

To obtain the programmed form of the equation of motion for 
the stator, the previous derivation is followed, but with 
U = 0, Hg becomes H, W becomes V, and B becomes a. 

B-3 THE AREA RESTRICTION FACTOR Z 



The condition at the outlet of a blade row with boundary 
layers on both sides of the flow channel is shown in Fig. B-1. 
The flow is considered to be turbulent within the boundary 
layer while, outside the layer, the velocity of the flow is 
the theoretical velocity. .Assuming a power-law velocity pro- 
file, the velocity may be written. 



Vth 




(B-45) 



The mass flow rate exiting the blade row can be expressed as 



nn = Pth Wh 



S' 












CoSo^cl 



+ 



0 (B-46) 

where and represent the ideal conditions for an 

isentropic expansion through the blade row to the discharge 
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pressure Pd, which is assumed to be constant across the blade 
spacing. The discharge angle of the flow leaving the blade 
row is closely approximated by the expression [Ref. 1] 



CB-47) 

Inserting Eq. (B- 47) into (B- 46) and reducing yields 



o<^ - cos' 



-I 



O- 



S' 



CO^o<^ 



m “ /^th Vth ^ 






Assuming a perfect gas 

P Tth 






T 



)q pru 



T-ro - -TthI 

TVo ■" CTtd -"Hth 



(B-48 ) 



D.efining 



y 



E = 



iri 

Pj \ <r 



TO 



(B-49) 



(B-50) 



Equation (B- 49) can be written 

P 



I - y E 



Substituting Eq. (B- 51) into (B- 45) yields 






m - PrH ^ 






(B-52) 
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Using the displacement thickness given by 






I - (i-Ve) 



n 



rr\ 






(B-53) 



the mass flow rate can be written as 



rv^ - ^TH ^ 






CU 



(B-54) 



The loss coefficient, expressed in terms of average kinetic 
energy lost is 









rv^ ^ Yju 

^ ’ (B-55) 

where E is the actual kinetic energy of the flow, given by 



s 1 

jur 



E = pTh^TH ( ^ ^ ^ y 



CB-56) 



Substituting Eq. (B-51) into (B-56) gives 






TH 



,2- 



(X 












CB-57) 
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The energy thickness is written as 






1 3m 



The loss coefficient can therefore be written as 

r 

I- 1 ir 



CB-58) 



f-'- 



I- 2 



11 

cu 



(B-59) 



The area restriction factor Z, is the fraction of the flow 
area through which the uniform theoretical velocity would 
produce the actual flow rate, thus 

CU 



Defining the energy parameter (a form factor) as 



H 















using Equations (B-59) and (B-61), Eq. (B-60) becomes 



z- 



- 1 






(B-60) 



(B-61) 



(B-62) 



where ^p is the profile loss coefficient, 
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B-4. THE ENERGY PARAMETER, H*** 

In Equations (B-53) and (B-58 ) the denominator of the 
integrand is expanded using the binomial theorem, so that 



/ arv,\-< 2^ Htr\ 3 

)-|+ye/] )(b Ji 4>y£/2 



(B 



The integral of Equation (B-58) is now written as 



s 






o 1 " Xe 



3r>\ Tirv^ 2. 7/>» 

+- Xg /] -i- ye +• 



V/ 3 ^ , 'I 

^ ^ Xe /l 4-... 

which, on integration becomes 

3m 






1 






+ 



Xe 



<=//}- -L ^ 









/U+I 

Therefore, Equation CB-58) becomes. 



(3 



■Xe . Xe 



7nh-H 



(B 









- I - 

II no+l 



Zrr\H 



yg 

Crri-tl 



(i-Xe) 



yg 

lrr\M 



y| 



CB 



-63 ) 
-64 ) 



-65) 

-I- 

-66) 
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which, can be written as 




In a similar manner. 



(B-67) 




I 



-h 

/V\-H 



Yg 

Srv\-H 






-I 



(B-68) 

Substituting Eq. (B-67) and Eq. (B-68) into Eq. (B-61) , the 
equation for H*** used in the computer program is obtained: 



H 






I ^ te. ^ ^ I 

y-l ^ Srn-H S'm-H ~/rr\H llfr\+l 



^ • Xe 

4- + -h ^ +* 

Xa'f 



(B-69) 
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FIGURE B-1: BOUNDARY LAYER EFFECTS AT THE EXIT 

OF A BLADE ROW 
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APPENDIX: C 



PROGRAM SEGMENTATION ON THE HP- 1000 

Segmentation allows large programs to be run on the 
HP- 1000. The program is divided by the programmer into a 
main program and several segments, which are stored on the 
disc. Each segment and the main program are then compiled 
and loaded. When the program is executed, the main program 
and its segments are called into memory individually, and 
only as they are needed for execution. In this manner, a 
program can run in a partition which is smaller than that 
program's total size. 

When the main program has performed all executable state- 
ments, the first segment is called into memory by an EXEC 
call. The system then loads that segment from the disc into 
a memory block following the end of the main program. The 
process is illustrated in Figure C-1. Note; the main program 
plus the largest segment may not together exceed 29 k. Once 
a segment is in memory it can call another segment. 

When executing, any segment can call any subroutine which 
is attached to the main program. It was this feature which 
allowed the present program to be run. All subroutines were 
placed within the main program. In fact, the main program 
consisted of nineteen subroutines and functions. A segment 
may not return to the main program. Communication of data 
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-■:== 




between the main program and the segments is accomplished 
through a common block. 

The four segments of the present program are "MAIN", 
"SHORT", "PART 2" and "PART 3". The manner in which control 
is passed from the main program to the first segment and from 
the first segment to the second is as follows: 

BLOCK DATA 



END 

PROGRAfi THESS 
DIMENSION INAM (3) 

DATA INAM /2HSH, 2H0R, 2HT / 



CALL EXEC (8, INAM) 

END 

PROGRAM SHORT (5) 

DIMENSION INAM (3) 

DIMENSION NAME (3) 

DATA INAM /2HSH, 2HOR, 2HT / 
DATA NAME /2HPA, 2HRT , 2H2 / 



CALL EXEC (8, NAME) 
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END 



PROGRAM PART 2 (5) 

DIMENSION NAME (3) 

DIMENSION NAMR (3) 

DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 



CALL EXEC C8, NAMR) 

END 

PROGRAM PART 3 (5) 

DIMENSION NAME (3) 

DIMENSION NAMR (3) 

DATA NAME /2HPA, 2HRT, 2H2 / 
DATA NAMR /2HPA, 2HRT, 2H3 / 



END 

The "(5)" after the program name indicates that it is a 
program segment. Note the manner in which the program name 
is put into a data statement using the Hollerith notation. 
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FIGURE C-1: PROGRAM SEGMENTATION- ILLUSTRATION OF THE J4AIN 

PROGRAI4 CALLING A SEGMENT INTO LOGICAL MEMORY 
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APPENDIX: D 



RUNNING THE COMPUTER PROGRA.M 

If the reader is unfamiliar with the HP-1000 Computer 
System, references [11] and [12] should be consulted before 
attempting to run the program. 

D-1. DATA INPUT 

Using the editor, input the following data into segment 
"SHORT'*. 

1. Turbine operating conditions: referring to Table 

A- III, type in appropriate data in lines 66 through 69 and 
74 through 78. 

2. Special input data/program control parameters: 
referring to Tables A- IV and A-V, type in appropriate data 
in lines 83 through 98. 

3. Turbine geometry: referring to Tables A-1 and 

A-II, type in data for stator and rotor in lines 103 through 
186. 

D-2 COMPILING THE PROGRAM 

1. To compile the main program type: 

:RU,FTN4,MAIN: : 25,- ,- 

2. Compile the first segment: 

:RU,FTN4, SHORT: :25,- ,- 
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3. 



Compile the second segment 



:RU,FTN4,PART2: : 25,- 
4. Compile the final segment: 
:RU,FTN4PART3: :25,- 

D-3. LOADING THE PROGRAM 
Type 

: RU, LOADR 
Tap return key 
Will display 
LOADR : 

Type 

OP, LB 

Will display 
LOADR : 

Type 

:RE,%MAIN: :25 
Will display 
LOADR: 

Type 

:RE,%SHORT: :25 
Will display 
LOADR: 

Type 

:RE,%PART2 : : 25 
Will display 
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LOADR: 



Type 

:RE,%PART3: :25 
Will display 
LOADR: 

Type 

:END 

After the end statement, the loader will display that the 
program is ready for execution. 

D-4 RUNNING THE PROGRAM 
Type 

: RUN, THESE 

The program will be executed and no further action by 
the operator is required. The computed pressure ratio of 
each iteration of the outer loop of the program is displayed 
on the screen as it is calculated. The operator therefore 
has some idea where in the iteration process the computer 
program is executing. 
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APPENDIX; E 



DISCREPANCIES IN MACCHI’S PROGRAM 

1. Main program, lines 21 and 22; the value of ICL has 
not yet been read. 

2. Main program, lines 163-166; the Traupel method of 
calculating gas outlet angles does not take the Mach number 
into consideration. However, in lines 163-164, the program 
is attempting to draw a parabola through points which repre- 
sent outlet angle as a function of Mach number. 

3. Main program line 281; the calling of subroutine 
SEINE is questionable. Parameters are transferred to that 
subroutine, but many of them have not yet been defined 
(HE, DHEDX , WPER2 , DSDXl) . These undefined variables will 

be set equal to zero by the IBM 360 and 370 computers. Thus, 
in line 10 of subroutine SEINE, the value of DWDX will be 
zero and in line 17, division by zero will occur and the 
execution of the program should cease. 

4. Subroutine ROTORI lines 22 and 26; the stator radii 
are used in the calculation whereas the rotor radii should 
be used. 

5. Subroutine ASOSI, line 107; the correct Fortran code 
is 

ZETAPS(I) = .5 * ZETAS(I) 
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6 . 



Subroutine AL0S2, line 121; the correct Fortran code 



is 

ZETAPRCI) = .5 * ZETAPR(I) 

7. Subroutine AL0S2, lines 123-126; the stator radii are 
used in the calculation whereas the rotor radii should be used. 

8. Subroutine ANGAIN, line 14; the correct Fortran code 
is 

AO = ATAN(1. -XCL/H*CH*C0SCANG1)/C0S(ANG2) * 

TAN(ANG2) + XCL/H*CL*COS (ANGl) /COS (ANG2) * 

TAN (ANGl) 

Note: Since reference [2] was published, Professor Macchi's 

program has been further developed by Professor Macchi under 
private sponsorship [Ref. 13] . The new code however, is not 
generally available. 
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APPENDIX F 



COMPUTER OUTPUT 
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APPENDIX: G 



COMPUTER LISTING 



240 



T=00n04 IS ON CR00025 USING 00147 BLKS R=0 000 



0001 
U002 
0 0 03 

0004 

0005 
00 06 
0007 
0 0 OB 
0009 
0 010 
0011 
0012 

0013 

0014 
OUiS 
0016 
0 01/ 
0018 
0019 
0030 
0021 
0 022 

0023 

0024 

0025 

0026 
0 0 27 
0028 

0029 

0030 
C031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 
0 039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 
0859 
0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 






c 

c 



BLOCK DATA 

C0NH0N/ABA/BAi7 , BLEX 
CCsMhON/C‘5R/CQSL( 1 0 ) 

COhMON/TOl./TOLl ,TnL2,T0L3,T0L4 
r.0nh0N/TR5/TRA5 

CnrlhON/r.AS/CP .GAM . FMMF. .ERRE . FXP L . EXP2 , 0 T 82 , OIS3 
r.nMM0N/cn7/icDR ,ino7, iiNC, lAi.icd, IAN, icon' 

COMLON/MAC/IN 
GilMMON/IW i/TND, 
rOMMON/AUS/XCL 
rOMMON/CSS/CJ .G,Oi 

CnrMON/UARl/RC( 10 > .RS0LD2 .RSnLD3.RSOLD4 ASFO.RSFO.RRFO.ARFO, 

:i^RR( ID) ,RRm,D2.RR0LD3,RR0LD4.Cy,CK .UAlCH)-) ,DALF(10 ) ,DBET(10 > ^ 

■^cASF AMS.Bi(20) 

c6mmON/C AR2/B6(20) , 7R , ZS . AR F . B2 ( 20 ) . PR , AMR , OUl ( i 0 ) 

COMhON/VAR3/P TE(lO) , RS t . RS3 . RS5 , T2 ( i 0 ) 

C0MM0N/UAR4/BR1 .BR2 .'RRS'.RRl '> R R3 , RR5 , VA2 ( 1 0 ) 

C0MM0N/VAR5/PRA1 i( tO ♦ .RlNCI (10 ) ,ALF Ali(iO) , BETA 1 1 ( 1 0 ) ,ZETA1 ( 1 0 ) , 
*L'2(iO) ,AI FA22(10) .HFTA22(10) 

CDMrtON/OAR6/PT2(iO ) > TT2( t 0 > , PTi ( 1 0 > , DEl.H( 1 0 ) , ALFA2 ( 1 0 ) ,0U2( 10 ) , 
tUR2( 10 ) , T2S(10 ) ,T2IS(iO) 

C0MM0N/UAR7/TTIS(5) ,BFTAT(5) ;ETAT(5) ,ETAl(iO) , ET AS ( 1 0 ) , ETAR ( 1 0 ) , 
KRSTARdO ) .AKXS(10),P8TR(i0) 

CUMM0N/VAR8/DR1 (iOK AMUKiO ) >AMU2(iO ) ,BFTET(10 ) ,PRAT1T(10 ) ,AMR2C 
:4C10) , YS( 10 ) XI ( 10 ) ,a 6FA1. ( i 0 ) , 7FTAPS ( i 0 ) ./JPFRl (10) . YR ( 1 0 ) , X2 ( 1 0 ) 
C0MM0N/UAR9/ZETAR < 10 ) ,ZETAPR (1 0 ) , AS ( 1 0 ) , AR ( 1 0 ) , SI 1 ( 1 0 ) , SI2 ( 1 0 > , 

( 10 > .DSDXl (10 ) 1( 10) ,HE( 10) 

COMMON/VARlO/UlJi (10 ) .DHEDXdO ) ,DSDX2(iO) , R I i d 0 ) , R 12 ( 10 ) , 
*813(10) .RI4dO) .RIdO) .SRldO),SR2dO) 



.AAdO ) .Sl^ dO) ■ PRAT2dO) ,UPER2dO) 
,.SSdO),ALFAdO) 



COMMON/UARii/YOLDd 0) 

*r'WDXdO) .TllSdO) .PPAT3dO) v r..., , 

C0Mn0N/VAR12/BEl Ad 0 ) ,DELR( 10) .WPtKO( 10) .ZFTASdO) ,7ETARi(20) , 
*Z£TAR3(20 ) , ZETAR5 ( 20 ) , R i ( 20 ) , Ai ( 20 ) , Tt 0 ( 20 ) 
r.CiMMQN/UARl3/STi (20 ) , IRR ( 20 ) , R? ( 20 ) . A2 1 20 ) , RING (20 ) ,DR(i0> , 

*BETO (10), ST AL 1 1 (1 0 ) , AREA2 ( ^ 0 ) , L‘R 1 ( 1 0 ) 
r:0hM0N/VAR14/ULBM,PRAr5,0MEG 

COMttON/ALi/ALFAl (10) , U1 ( 1 0 ) . TTO , RPM , R S (1 0 ) , SI , TNI , H , D , Cl , T1 d 0 ) ^ 
•YP 1 ( 1 0 ) , TO , TE I . ALI , BFSP , XX , AfQG20 , AMS 1 ( 1 0 > , S , TN , C , Tc , AL , SO , TNO , 

*r.O TEO,n(iO) ,Dll,DlO,ni*>l .D20,ANr;2I ,ALFAX,T1 J .PTO . alo . AM(i 
C0MMaN>AL2/BeTA2d 0 ).,HETAl (1.0 ) , BET AO( 1.0) :W2( 10) , TtR (1 6 ) , U2 d 0 ) > 
*SXR,1NIR .HR ,DZ,C1R ,TXPC,SZ,TNR .r.R .SOR , TNOR , COR . Al_T R , ALR , ALOR . 
*P2.(iO ) ,1402(10 ) ,W1( (O) ,TEIR, ter; TEOR,Di.TR, DIOR, BETA7,BETAI^ANH^ 
*TIR-.IR ,TDR ,STALId0) 

COMHON/ARF/REE 

COMMON/TRA/XPOl (5,8) ,XP02(6,8) , AI,F 1 ( 8 ) , Al. FOl ( 5 ) ,AI..F02(6) , 

*YdO) , Yt (10 ) ,0(6) ,RX(30 ) ,RY(30) ,IR(30 ) ,Z(6) ,Ct <4,8),C2(4,8) 

DATA ALFI/ i6 . , 15. ,20 . ,25 . ,30 . ,40 . ,50 . ,60 ,/ 

DATA ALFOl/ 40 , , 50 ! , 60 ! , 70 . , 80 ' / 

DATA ALF02/ 80 . , 90 . , 1 0 U . , 120 . , 150 . , 170 , / 

DATA XPOl/ . 0570 ,. 0465, . 0440, . 0428, . 0424 



DATA XP02/ 



. 0530 , 
. 0495; 
.0475, 
. 0440 , 
.0420 , 
, 0420, 
. 0420 , 
. 0424, 
. 0323 



041S, . 
. 0380; 

. 0355, . 
. 0335, 
.0312, . 
. 0300 , 
.0300, . 
, 0422, 
0320 



0350 , 
0312, 
029S, 
0273, 
0224, 
0213; 
0206; 
0420 . 
0318 



, 0330 , 
. 0296; 
, 0267, 
0245, 
. 0205, 
, 0i8i; 
,0155, 
0402, 
0295 



. 0323, 
. 0285, 
. 0250 , 
. 0225, 
.0183, 
.0152", 
. 0125/ 
0313 
0200' 



0000 , 

0000 . 



;0283; 0280", 027*5 ; ! 0250 0 3.43 ioOOO", 



. 0250 
, 0225', 
,0183, 
, 0150 , 
.0125, 



0246, . 
0216; 
0170 , 
n 1 36 ; 
009?; 



0242 , 
0203, 
0154, 
0104; 
0 0 73, 



0208 . . 0070 , 

0168. . 0000 : 

01 06, . 01 no , 
0050 ntsd 
0 00 0 . n?n , 



0000 
0000 
0000 , 
0000 



. . ..w ^ ..r.. , .000 0 / 

DATA UIS2, 9IS3, CP ; FMMF , gam/. 0^10013, .00001^^, .24,28.97,1.4/ 
DATA G, CJ, EXP 1, EXP?, ERft E/32. 174,778 . 1 6 , 3 . 5 , . 2857 , 53 . 3459/ 
END ^ ^ > > 

PROGRAM THESS 



DIMENSION INAM(3) 
COmMON/ABh/BAI 7 , BLEX 
COMMON / CUR /COSLdO) 

COMMON/ rnL/TfJLi , T0L2 , T0L3 , T0L4 
CnrtM0N/TPS/TRA5 
CnMMitN/i:A8/r.P .GAM ,FMMF ,LRPF 



EXP,? ,UT82 . VTS3 



?41 



fif)7Q 

0080 

0 081 

0082 

0083 

0034 

0088 

0036 

0087 

0 0 88 

0089 

0090 

0091 

0092 

0093 

0094 

0098 

0096 

0 097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0103 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

U132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

10157 

0158 



CnMM0N/C07/IC0R , ICOZ . I INC , I A I , I CL , IAN , I CON 
- rXirlhON/hAC/XN 
COnhON/IWI/TND, INZ> lUR 
CO.-nhON/AUS/XCL 
C0r,M0N/CS5/C J . G , a 1 

COhhON/UAR 1/RC( 1 0 ) >RSni.D2 ,RSni.D3 ,RS0IX)4 .. ASFO .RSFO . RRFQ, ARFO, 
,RR0I.D2.RRrjl.D3.KR0LD4,C7,CK A'AKIO) .DALFdOl ,DEET(10) , 

ItASr ,ArtS,Bl(20^ 

CO(nhON/VJAR2/H6(20 ) ,7R 2S • ARF 32 ( 20 ) . PR , AMR , V^Ul ( 1 0 ) 
COMMON/yARX/PTEOrj) ,R5l. . RS3 , P55 , T2 ( i 0 ) 
rX‘Mhi:iN/L/AR4/HRi .&R2 i PR-Tl'^RRl : RR3 . RR5 > VA2 < 1 0 ) 

COmMON/VJARS/RRAI 1 ( nv. .RINCI (1 0) ,ALf- All (10) .BET A1 1( in ) ,ZETAi(10) , 
)kl»2(10),AI.FA22(i0> RFTA22(10) ^ ^ ^ 

COMnQN/L^AR6/PT2(l0) ,TT2(10 ) ,PTi( 10) ,DF!_H(10) , ALFA? ( 1 0 ) , V^ll2 ( 1 0 ) , 
'}cWR2(10) ,T2S(iU) ,T2T<o(10) 

n0MM0fVUAP7/TTlS(5> .PFTAT(5) ,ETAT(5) ,FTAI(10) ,ETAS(10) ,ETAR(10) , 
:*^RSrAR(i0),AKIS(10) ,PSTR(10) 

COnMON/yAR8/DRl(103 ,AMU1(10) ,Ahy2(l 0) , BFTFT ( 1 0 ) , PR ATIT ( 1 0 ) , 
:*<AMR2(10) , YS ( 1 0 ) , X t (1 0 ) ,AREA1( 10) ,ZETAP8(1 0) , UPEP 1 < 1 0 ) , YR ( 1 0 ) , 
XCX2C10) 

COMMC)N/yAR9/ZETAR( 1 0) .ZtTAPRdO) , AS ( 1 0 ) , AR ( 1 0 ) , SI 1 ( 1 0 ) , SI2 ( 1 0 ) , 
#Si (10) ,D5T>Xi(l0) ,un ( to) ,HE(10) 

COMMON/‘JAR10/UU1 (10) .DHEDXd 0 ) ,DSDX2( 1 0),RI1(10),RI2(10), 
^RI3(10);RJ4(lO),Rr(iO)>SRl(in).SR2(10) 
rnMhON/yARM/YOLDdO) .AA(10) XR(10) ,PRAT2(i0) ,UPER2(in) , 

.4^J)1‘IDX( 10) ,1113(10) ,PRAT3(10) ,33(10) ,ALFA(iO) 
r.rjMMClN/yAR12/BE1 A( j 0 ) ,DELR (10 ) UPERQ( 10) ,ZETA3( i 0 ) ,2ETARt(20) , 
)KZhTAR3(20) ZErAR5(.?0) ,Ri(20) .Ai (20) ,T10(^!0) 

C0hM0N/yAR13/&Tl (20) J RR(2.0) ;R2(20) ,A2(20) ,R 1NC( 20 ) ,DR ( 10 ) , 
YPFrOdO) ,3TALII(l0) ,ARFA2(in),yRl(lO) 
CnMM0N/yAR14/WLBM>PRAT5,0MEG 

COMMON/ALj /AI.FAK j 0) A'lUO ) .TTO ,RPM ..R3( 10) ,SI ,TNT . H , D , Cl , T 1 ( 1 0 ) , 
:ic?l (10) ,TO,TEI ,ALI , PF3P , XX , hNC? 0 : AMSl ( 1 0 ) , 3 ' TN C , TF , AL , SO ; TNO , 
i:r.O ,TEO,li( 1 0) ,6lJ. ,t>20 , ANG21 . ALFAX , 1 t T , P tO AiX.AMC 

C0MM0N/AL2/B^TA2( 10 ) ' BETAl ( lO ) ,PETA0( 10 ) ^J2^ if^ ) . TTf.:( lO ^ ,U2( 10 ) , 
*3 Tr inir .HP ,dz,cir >ttpc,sz,tnr ,cr ,sqr ,tmor ,cor ,Altr , ALR , ALOR . 

Y:P2(i0 ) ,iJn2(i0) ,W1 ( I 0 ) , TEl R , TER , TEOR , Di TR . DIOR , BETAZ BETAT , ANM , 
)tcTIR,TR,TOP,STALl(10) 

COMMON/APF/REE 

COMMON/TRA/XPOl (5,8) >XP02(6,B) »A1-F1 (8) , ALF0K5) ,ALFn2(6) , 
«Y(10)vYl(l0),O(6),RX(30),RY(30),IR(30),Z(6),Cl(4,8),C2(4'8) 

DATA fNAM /2HSH , 2H0R , 2HT / 

TRA5=1 

XX=1.25 

CArUL EXEC(8,INAM). 

END 



SUBROUTINE CHAN O TO . AMC , P TO , RC , UlLBM , UCHAN , WPERO ) 
DIMENSION RC(IO) ,UPeRO(iO) 

COMMON/GAS/CP , GAM, EMME, ERR F, EXP 1 ,EXP2, 0132, yiS3 

C0MM0N/CS3/CJ,G,0i 

TC=TTO/< 1 .-KGAM-i . ) /2 . *AMC*AMC ) 

yC=SORT(GAh)*cERRE*G*TC)*AMC 

PC-pTO/d , + (GAM-l , )*aMC*^1c2/2. ) *:<cEXP i 

RHO=PC/ERRE/TC 

AREA=3 . 1 41 6*(RC(5)**2-RC( 1 )**2) 

ULBM=RHOxAREA*VC 

wnHAN=WLBM/(PTOi*cSQRT(G/ERRE/TTO) ) 

UPEROd ) = 0 . 

UPER0<2)=.25 

UIPER0(3) = ,5 

UPER0(4)=,75 

UPER0(5)=1 . 0 

RETURN 

END 



SUBROUTINE STATR ( ALF A1 , X , TTO , PTO , AM , T , P ,Vi ,UA1 , SI 1 , SI2 , Y , S , DSDX 
)*:yUl,PRAT,TlIS,SS,DALF,RSF,DELR,CL,CK,7FtAP3,R,RSl ,RS3,RSS, 
*7FTA,DR ,7ETAS,AMS ,NS,URi) 

DIMENSION ALFAl d 0 > , X ( 1 0 ) , T d 0 ) , P ( 1 0 ) , U 1 d 0 ) ,UA1 ( 1 0 ) ,SI1 d 0 ) , 
•1cS12(l0) ,V(10) ,DSOX(^ 0) ,VLli(in) . P R AT (1 0 > . Tl 1 3 d 0 ) , S3 ( i 0 ) , 3 ( 1 0 ) , 
^DALDXClfO ,AI.PAd0) ,AL.FAMd0) ,DALF(10) ,AM3(10) , DAIFDX ( 10 > , DFLR ( 

.^1 0 ) ,ZETA3d0 ) ,ETA( 1 0) ,2ETAPSdO ) , R (10) . 7FTA( 10 ) , DR ( 1 0 ) , UR 1 d 0 ) 
C0MM0N/GA3/CP,GAf1 ,EhME,ERRE,EXPi,EXP2,0XS2,UIS3 
rOHhON/C:P3/CJ^ G,Oi 
nLtMMON/IUI/IND,INZ,IUR 
C8=0 . 0 



C9=0 . 0 

01=2 . :4cC.T*G4fCP 
n2=UAl ( 3 > ^4C2 / ( Ql^TTO ) 



> 



2^2 



0159 

0i60 

OiSi 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

Ol77 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0136 

01B7 

0188 

0.189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

020R 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 



. DO 303 1=1,5 

TF(f<<I)-R83) 300,301,302 

30 0 ZFTHS(I)=7ErA(l)4.( (R( X)-RS1 )/(R83-RS.1 ) > ( 7FT A (3 ) -ZFTA ( 1 ) ) 

AI.FAl ( I )=A!.FA1 (3) + ( (R ( I )-RS3)/(RSl-RS3) )xDALF( 1 ) 

7FTAPS<T)=ZETAPS(l) + ( ( K ( I ) -R81) / < RS3-RS 1 ) ) * ( ZF.7 APS ( 3 ) «ZET APS ( 1 ) ) 

00 TO 303 

301 ZFTAS(I)=7ETA<3) 

ALFA1(I)=ALFA1(3) 

00 TO 303 

302 ZETAS( 1 )=7FTA(3) + ( ( R ( I) -R S3 ) / < RS5-R S3 ) ) :^ ( 7ET A ( 5 > -ZF.TA ( 3 ) > 

ALFA! ( I )=ALFAl<3) + < <R(I)-RS3)/(RS5-RS3) ).icOALF(5) 
7ETAPS(I)=ZETAFS(3) + < ( R ( 1 ) -RS3 ) / ( R35-RS3 ) ) ( 2ETAPS < 5 ) -ZET APS ( 3 ) ) 

303 CONTINUE 

DO 305 1=1,5 
ETA(I)=1.-ZETAS(I) 

M=I-1 
N=I + 1 

IF(l-l) 307,307,309 

307 DALFDX(l)=(ALFAi(2)-ALFAl(l) ) / ( X ( 2 > -X < 1 ) ) 

GO TO 315 

309 TF(I-5) 311,313,313 

311 DALPDX(I)= .5*(< ALFA1<N)-ALFA1(I> )/(X(N)-X<I) )+ ( ALFAl ( I ) -ALFAl ( M) )/ 
*(X<I)-X(M) )> 

GO TO 315 

313 DALFDX(5)=<ALFA1(5)-ALFA1(4) )/(X(5)-X(4) ) 

315 TAN1=-2.*TAN(ALFA1(I)) 

PROD=TANi*DALPDX( T ) 

SINS0=-2 . i^cSTN( ALFAl < I ) )**2/X( I ) 

SH(I)=PRnD+SINSQ 

305 CONTINUE 

304 DO 332 J=l,5 
IF(J-l) 306,306,310 

306 IF(NS-l) 317,310,310 

317 DO 308 1=1,5 
SS(I)=0. 

308 SI2(I)=SI1(I) 

GO TO 318 

310 DO -312 1 = 1,5 

AA=C2* Y ( I) **2/C0S ( ALFAl ( I ) ) 

AB=(1 .-AA)/< 1 .-AA/FTA(I) ) 

312 S(I)=ALOG(AB) 

314 DSDX(1)=(S<2)-S(1 ) ) / ( X ( 2) -X ( 1 ) ) 

DSDX(2)=0 .54C(DSDX(1 ) + ( S <3 ) -S ( 2 ) ) / (X ( 3 ) -X < 2) ) ) 

DSDX(3) = .5)*c( (S<4)-S(3))/(X<4)-X(3))-^'(S(3)-S(2))/(X(3)-X(2))) 
DSDX(4) = .5:*C( (S<5)-S<4) )/(X(5)-X(4> > + (S(4)-S(3) )/<X(4)-X(3) ) ) 
DSDX(5)=(S(S)-S(4) )/(X(5)-X(4) ) 

DO 316 1=1,5 
IF(NS-l) 319,321,321 

319 SS(I) = ( 1 .-COS(ALFAl( I) ):|t)|c2/ (C2*YC 1 )*:<c2) )*DSDX( I ) 
no TO 316 

321 SS(1) = ( C-COS(ALFAl ( I) ) :^c*2/ ( C2«Y ( I ) )«{*2 ) ) +STN ( ALF A1 ( I) ) *)«(2+C0S ( AL 
:tcFAl ( I ) )*)<(?)lc(CL**2+(DR ( I)/2. 0 ) *)T2 ) > *DSDX ( I ) +COS ( ALFA! < I ) 

%CK%2 . :<CRSF*DELR < I )/CL«*2 

316 SI2(I)=SS(I)+SI1(I) 

318 SUM1=(SI?(1)+S12(2) ) ( X(2 ) -X ( 1 ) > /4 , 

SUM2=(ST2(2)+SI2<3) )«(X(3)-X(2) )/4 , 

SUM3=(SI2(3)+SI2(4) ))K(X(4)-X(3) >/4. 

SUM4=(ST2(4)+SI2(5) ))K(X(5)-X(4) )/4 . 

EN2=-SUH2 

FNi=-SUM2-SUMl 

FN3=SUM3 

FN4=SUM3+SUH4 

Y(2)=EXP(EN2) 

Y(i )=EXP(ENi) “ 

Y(4)=EXP(EN3) 

Y(3)=l . 0 
Y(5)=EXP(EN4) 

IF(IND-l) 332,323,323 

323 IF(J-l) 324,324,320 

320 IF(J-3) 322,324,322 
32? JF(J-S) 332,324,332 

324 UIRITE(6,326) 

326 FORhAT (/57H SLINE C8 C9 ITERATION I ^ ALFA I^DSDX , I'TOTAL 
^ Y ) 

"nn 330 1=1,5 

328 FORhAT ( 14 , F4 . 2 ,F4 . 2 , 19 , F 12 , 4 , F9 . 5 , F 9 . 4 , ?F8 . 4 ) 

330 WRnE(6>328) 1 , C8 , C9 , J , Si 1 ( I ) , SS ( 1 ) , SI2 ( I ) , Y ( 1 ) , AI.FAl ( I > 

332 CONTINUE 

DO 334 T=i ,5 

•JAi. ( I)=UA1(3).1 cY( I ) 



243 



0239 

024U 

0241 

0242 

0243 

0244 

0245 

0246 

024 / 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

025/ 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

027/ 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 



334 

336 

352 

356 



VUl(I)=VAi ( t ))*:TAN(A)..FA1(I ) ) 

U1. ( I)=^A1 ( D/CnSCAl.FAl (I ) ) 

URi(1)=-Vr1 (I):<cDR(I)/2./CL 

'/I < 1) =5QR T ( U 1( 1 > 1 ( I ) +VR t ( T ) WR 1 < I ) ) 

T( I )=TT0-U1 ( 

riisa )=TTn-nio-T( T > )/eta(i > 

P < 1 )=PTO*( T1 CS( I )/Tin))|c.¥EXPl 

PRAKI)=P<I)/PTO 

DO 352 1=1,5 

AMR ( I ) 1 < I ) /SQRT < GAM>KERR < I ) ) 

CONTINUE 

RETURN 

END 
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SUBROUTINF ROTOl <UU1 . VAi . RPM ,U . BET A1 , HF , TTE ,PTF , X2 . P 1 .T1,W1 ,U 
«X1 ,RS , ZFTAR ,ZETAPR . RR . DHEDX . DSDX , S Oi.2 , OMFG , BR i ; BR2 , BR3 , FSl ;FS2 > 
)«7FTA,&6,RS^ ;RS3 .RS^ .BF.TO,BT(^^I.I ,BINC,UR1 ) 

DXhENSlON ,UA1 (5) ,U(5) . BET A1 ( S ) . HF ( 5 ) > TTF < 5 ) ,PTE(5) , 

:icX2(5) ,P\(5) , T1 < 5 } , U1 ( 5 ) ; WUl ( 5 ; > XI ( 5 ) .R8(5) ,ZETAR (SJ , 

>K7FTAPR (5) ,RR<5) .DHFDX<5> . DSDX < 5 ) > S < 5 ) . U2 ( 5 ) ,2ETA<5) , 

.icVRi<5) .B6(10) ,BtTn<5) ,STALI<5) ,BINC<5) 

COMMON/CSS/GJ >G , 01 

COnnON/GAS/CP .GA*i .EMME,ERRE,EXP1 ,EXP2>UIS2 ,UIS3 
COMnUN/IUI/TND, IN7, lUR 
C=2 . )K3? . 1 74*778 . 16*CP 
0NEG=RPh*3. 1416/30 . 

DO 520 T=i,5 
U(I)=0MEG*RS(I)/12. 

U2(I)=U< I )*RR(1)/RS<I) 

UUl (I)=UU1.(T)*-U(I) 

PETA1<I)=ATAN(UU1 (T)/UA1(I) ) 

W1 < I)=UA1 <.T)/CCS(BETA1< I) ) 

U 1 < I ) =SQR T < UR 1 < 1 ) *UR 1 ( I ) +W 1 ( T ) *U1 < I ) ) 

TTF(I )=T1 ( X)+Wl ( I )*Jk 2/C+<U2( T )**2-U<I)**2)/C 
PTF<I)=PiCI)*(TTE(I)/Tl<l) )**EXP1 
HE(I)=TTE<I)*.24 
IF(RS(1)-RS3) 512,514,516 

ZETAR<I)=ZETA(1 )+(RS<I)-RSl)/<RS3-RSl)*(ZETA(3)-7ETA(i) > 

GO TO 518 
ZETAR<I)=ZETA<3) 

GO TO 518 

ZF.TAR(I)=7ETA(3) + <RS(I)-RS3)/<RS5-RS3)*(ZETA(5)-ZETA<3) ) 

ZETAP R < I ) =ZETAR < I ) /2 , 0 
CONTINUE 

DSl)X(l) = <S<2)-S(i) )/<X2<2)-X2<i) ) 

DSDX(2)=0 .5*<DSDX( 1 )4-<S<3)-5(2) )/<X2<3)-X2<2) ) ) 

DSDX(3)=0 . 5*< (S(4)-S(3) )/(X2<4)-X2<3) )+(8(3)-S(2) )/<X2(3)-X2<2) ) ) 
D8DX('5) = <S(5)-S(4) )/(X2(5)-X2<4) ) 

PSDX<4)=0 . 5*<DSDXt5)+(S<4)-S(3) )/(X2< 4)-X2(3) ) ) 

DHEDX( 1 ) = <HE(2)-HE( 1 ) )/ ( X2 < )-X2( 1 ) ) 

DHEDX(2) = . 5*(DHFDX( 1 )-^(HE(o)-HE<2) ) /<X2(3)-X2(2) ) ) 

DHEDX(3)=0 ,5*< (HE(3)-HE(2) ) / ( X2 < 3 ) -X2 ( 2 ) ) + <HE<4 )-HE< 3) ) / 

* (X2<4)-X2<3))) 

DHEDX<5)=(HE<5)-HE(4) )/(X2(5)-X2<4) ) 

DHEDX(4)=0 .5*<DH£DX<5)+<HE(4)-HE<3) )/<X2(4)-X2(3)) ) 

522 CONTINUE 
RETURN 
END 

SUBROUTINE R0T02 < BFTA2 , HE , DHEDX , DSDXI , DSDX2 , UA2 , WU? ^W2 , UU2 , U2 » 

*X2,U,YR .ZETAR.RIl ,RI2,RI3 ,RT/kRT,<^R 1 .SR2.AA,Sf^,TTI^,PT£,T2,P2,PRAT2 
*,T2S, INDS.DBET ,RRF , DFLR ,CL > CK . DR ,RRi . RP3 , RP5 ,NB , UR2) 
f>IMENSIOW BETA2<5) !hE<5) .DHFDX<5) ,DSDXU 5) DSDX2 ( 5) , UA2(5) , 

XUU2(5) ,W2(5) ,UU2<5^ , UP < 5 ^ , X? < 5 ) , U( 5 ) , YR < 5 ) , ZET AR ( 5 ) , 

*RI1<5) ^RT2(5) ,RI3(5rRT4<5> ,RI(5) . SR2<5) , YOLDCS) , 

*AA(5) ,SR (5) ,TTE<5) ,PtF(5) . TP < 5 ) . P2 ( 5 ) , PR ATP (5) ,TPS(5) , 

*DBET(^) ;BFTAf1(5>,AAR (5) ,DBETDX (5 ) ,B£TA <5) , DELR (5) ,RI5<5) , 

*DR(5) .R<5) vUR2<5) 

COMhON/TOL/TOLl >TGL2,TOL3>10L4 
CnNMON/CSS/CJ,G.Ql 

CnOrtUN/GAS/CP .GAh,EMMF,ERRE>EXPl ,EXP2,UIS2,UIS3 

CONMON/IWI/IND, INZ, lUR 

INDS=0 

TNDS1=0 

C=2.*G*CJ 

Q1=C/UAP(3)**2 

DO 274 T=1 ,5 

IF<k(I)-RR3) 270,?71,P73 

270 BFTA2< I )=BETA2(35 + (R ^ X ) -RR3 ) / <^ RR 1-RR3 ) *DPET (1 ) . 
r.n TO 274 
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271 8ETA2<1)=BETA2<3) 

GO TO 274 

273 BETA2( I )=BeTA2(3)-KR ( I) -R R3 ) / ( RR 5-RR3 ) < 5 ) 

274 CONTINUE 

BBETDXC 1 )=(BETA2(2)-BETA2( 1) )/(X2(2)-X2( i ) ) 
DHETDX(5)--:-(BETA2(5)-BETA2(4) ) / ( X2 ( 5 ) -X2< 4 ) ) 

DO 280 I=2>4 
M=T-1 

280 DBETDX(I)=,5:«< (BETA2(N)-BETA2< I) )/(X2<N)-X2(I))+(BETA2(I)-BETA2 

>/(X2( I)-X2(M) ) ) 

DO 10 1-1,5 

200 TANl=-2.:icTANCBETA2<I) ) 

PROD=TANl5f^DBETDX( T) 

5INl=-2 . *8TN < BE! A2 ( I) ) JK*2/X2< I ) 

R 1 1. < I) =PROD+SINi + DSDXl ( I ) 

SRl(I)=-4 . ^U(3))*tCnS<BFTA2(I ) )^STN(BETA2( X) )/<UA2(3)*YR( I) ) 
SR2(I)=2.)KU<3)*U( I ).YCnS(BETA2( I) ) )|{»c2/ < UA2 ( 3 > ( I) *Xc2 ) 

YnLD(l>=YR(I) 

AA(l)=(UA2(3)*YR<T)/CnS(BETA2<I) ) )*t2/(n*HE(I) ) 
RI3(I)=<C^C0S(BETA2(I) )*)C2/( ( UA2 ( 3 ) *YR ( I > )*)«c2) )*DHEDX( I) 

TF (INDSl-i) 10>250,250 
10 CONTINUE 

281 IF(lND-i) 201,282,282 

282 URIT£(6,121) (RIKI) , T = l,5) 

121 F0RMAT(/23H CONSTANT INTEGRAND l-5> 5F8.5) 

WRITE (6,122) 

122 FORhAT(/60H SEINE INDSl GRAD S INT2 INT3 INT4 INT 

UAL) 

201 DO 20 J=i,13 
no 30 X=t'5 

AA(1)=AA( I )*<YR( I )/YOLD(I ) )2H*2 
ANUM=i.-AA(I) 

' ADEN=i.-AA(I)/(i.-ZETAR(I)) 

AB=ANUM/ADEN 
TF (AB) 130,130,30 
130 TNDS=i 

GO TO 150 

30 SR(I)=ALnG(ANUM/ADEN) 

DSDX2(1 )=(SR(2)-SR(1) ) /(X2( 2 ) -X2 ( 1 ) ) 

DSDX2(2)=0 . 5)lc<DSDX2( 1 ) + (SR(3)-SR (2) ) / ( X2 ( 3 ) -X2 ( 2 ) ) ) 

DSDX2(3)=0 . 5*( (SR(3)-SR(2) ) /<X2( 3 )-X2(2) )-»•( SR (4) -SR (3) )/ 

% (X2(4)-X2(3) ) ) 

DSDX2(5)=(SR(5)-SR<4) )/(X2(5)-X2(4) ) 

DSDX2(4)=0 ,5»(DSDX2(5)+<SR(4)-SR(3) )/(X2( 4)-X2(3) ) ) 

DO 40 1=1,5 

SR 1 ( 1 ) =SR i ( I ) »YOLD ( I ) / YR ( I ) 

SR2( I)=SR2(I )*(YOLD( I)/YR( I) )**2 
RI2(I)=SR1( T)-SR2(I) 

RI3( I )=RI3( T)*(YOLD(I )/YR(X) 

IF(NS-i) 31,32,32 

31 RI4(I)=DSDX2( I)-(DSDX1‘( I)-hDSDX 2( I) )*Ci*HE<I) 

x(C0S(BFTA2(I) )/YR(I) )**2 
GO 10 40 

32 RI4(l)=-(DSDXi( 1)+DSDX2(I ) )KC1 :«HE( I ) * (CnS(BETA2( I ) )/YR( I) 

RI5(I ) = (DSI)Xi(I)+DSDX2(I) ) ( SIN ( BETA2 ( T ) ) 5K)«c2+COS < BETA2 < I ) 

**( (CL**?-»-(DR (I)/2. 0>**2)/CL5K)»:2) ) -COS ( BETA2( 1 ) ) ( 2 . *CK*RRF* 

*DELR(I) )/CL**2 
RI4(I)=RI4<n4.RI5(I) 

40 RT(I)=RI1 (I)+RI2( I)+RI3(I)+RI4(I) 

SUhl=(RI(l)-»-RI(2) )*(X2(2)-X2( 1) )/ 4 . 

5Un2=(RI(2)-».RI(3) ))lc(X2(3)-X2(2) )/4. 

Sim3=(RT(3)+Rl(4) ) X(X?(4)-X2(3) )/ 4 . 

SUN4=(RTc4)+RI(5) )*(X2(5)-X2(4) ) /4 , 

RN1=-(SI)M2-SUM1 ) 

EN2=-SUM2 
EN4=SUh3 
EN5=SUr13+SUM4 
DO 50 1=1,5 
SO Y0LD(I)=YR( I) 

YR(1 )=EXP(ENi) 

YR(2)=eXP (EN2) 

YR(3)=1 . 0 
YR(4)=EXP(EN4) 

YR(5)=EXP(EN5) 

NCOUNT=0 
DO 1001 T=i ,5 
XFCYRd ) .GT.2,0) YR(T)=2.0 
XF(YR(I) .1 T.0.2) YR(I)=0.2 
inni noNTTNUF 
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DO 110 1=1.5 

, TFST = A8S(Yf^l D(1)-YR<D) 

IF ( TEST-TP.L4 ) 110, 15. 0,119 
110 NCnUNT=NCni)NT + l 

IF(NC0Llf'^T-5) 119,140,119 

119 IF(IND-i) 20,120,120 

120 TF(J-3) 80,100,30 
80 TF(J-6) 90,100,90 
90 IF(J-9) 160,100,160 

160 TF<J-12) 20,100,20 
100 DO 60 I=i>5 

123 FORMAT ( 14 , 17 ,F10 . 5 ,5F8 . 4 ) 

60 MRITE (6,123) I , INDBl , DSDX? ( I ) , R I2( I ) , R 13 ( I ) , R J 4 ( I ) , R I ( I ) , YR ( I ) 
20 CONTINUE / > > ^ / 

140 DO 70 1=1,5 

UA2< I )=YR <I ))i(UA2<3) 

W2<I)=gA?(I )/L;0S(BETA2(I) ) 

UR2 ( I > =-UA2 ( I ) *DR ( I ) /2 . /CL 

U2 ( I ) =SQP T < U2 ( 1 ) ( X ) +UR 2 < I ) ( I ) ) 

70 T2<I)=TTF.<T)-W2(I)>»c)K2/(0.24:<cC) 

IFdNDSi-?) 251,149,149 
251 INDS1=INDS1+1 
DO 250 1=1,5 

AMR ( 1 )=U2( I )/SQRT (OAM*ERRE)»{GJ<(T2( I) ) 

250 CONTINUE 
549 DO 190 T=l,5 

WU2(I)=UA2d)«TAN(BETA2(I) ) 

VU2(I)=U1IJ2(I)+U<I) 

T2S(1)=TTE( I)-(TTE(I)-T2< I ) ) / ( 1 . -2ETAR < I ) ) 

P2(I)=PTE(I )*(T2S( I )/TTE(I ) ))!?*( GAM/ (GAM- 1 . ) ) 

PRAT2 0 )=P2(1)/PTE(I) 

190 CONTINUE 
150 RETURN 
END 

SUBROUTINE FLOUR ( PR AT . ZETAP .X , UI , PTE , P TO ,TTE , TTO , AS ,ZS ,RS , AR , ZR 
*RR,M,WChAN.UA,UPER,CnDE.Wt.BM,B,R .TIPC.A) 

DIMENSION PRAT (10) ,7ETAF(10) ,X(l6) ,UI(10) , PTE ( 1 0 ) , TTE ( 1 0 ) , 
*UA(10) .UPFR(IO) ,B(20) ,A(10) ,R(10) 

COMMON/CUR /COSLC 10 ) 

COMMON/TOi./TOLi , TOI-2 , TOL3, T0L4 
COMMON/MAC/TN 

COMMON/GAS/r.P ,GAM,EMME,ERRE,EXP1,EXP2,UIS2,UIS3 

COMMON/CSS/CJ,G,Qi 

COMMON/ APA/BAt7,BLEX 

COMMON/IUI/IND, IN2, lUR 

IN=20 

C=BLEX 



A(3) = B(l)+B<2)*R(3)+B<3)5KR(3))lc5tc2+B(4)5KR(3)JK)|(3+B(5)5lcR(3)#^4 
Fi=l ./(C+1 . ) 

F2=1./<3,*C+1. ) 

.F3=l ./<5.*C+i . ) ‘ 

F4=l./(7.)lcC+l.) 

F5 = l ./(9,)(cC + l . ) , 

F6=l./(11 .3KC + 1 . ) 

• PRATCR=(? ./(GAM+l , ) )**( GAM/ ( GAM-1 . ) ) 

PHICR = (2 . /(GAM+l . ) ):«:*c(i ./(GAM-1 . ) )*S0RT(2 . >»cGAM/ ( GAM+l . ) ) 

DO 420 1=1,5 

IF (PRATCR-PRAT(D) 400 ,402.402 
400 XE=1 .-PRAT(I)**( (GAM-1. )/GAM) 

GO TO 404 

402 XF = 1 ,-PRATCR)K»( (GAM-1 . )/GAM> 

404 XF2=XE**2 
XF3=XE)*c3^:3 
XE4=XE**4 
XF.INU=1 ./(XE-1 . ) 

HNUM=XFJNU+F2+XE*F3+XF2*F4+XF3t.F5+XE4)KF6 
H0EN=XEINU+Fl + XE4cF2+XE2*F3+XE3*F4+XE4)lcF5 

hstar=hnum/hden 

XI = (HSTAR-1 . )/<HSTAR-1 .+ZETAP( T) ) 

IF (PRATCR-PRAT(I) ) 406,408,408 
4 06 PHI=SQRT(2 . JUGAM/ ( GAM-1 . ) * ( PR AT ( I ) ** ( 2 . /GAM ) -P R AT ( I ) 5K* 

*< (GAM+i . )/GAM) ) ) 
no TO 410 
408 PHT=PHICR 

410 A(1)=B(5 )+B(2)*R( I)+B(3)*R(T )>fc)K2+B<4)5*cR( I)*)|C3+B(5)^R(I)#*4 
ARAT=A( T)/A(3) 

TF(M-2) 415,412,415 
412 TF(I-5> 41.5,414,414 

4 aRA!=ARAT+2% *3. 1 41 6*R < 5))<cTXPC/ ( 7R:*AR)tcRR.*: ( X ( 5 ) -X ( 4 ) ) ) 
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TF(IND-i) 120,416,416 
URT.T£<6 .41 8) XI.PHI.ARAT 

FORMAT (/SH XH=> F6 . 4 > 6H PHT=, F7 . 5 , 7H ARAT=.F6.4) 
Wj: ( I) = (PTF(I )/PT0)/5nRT(TTF( T >/TTO)*ARAT.'r.Xl4cPHI*C:nSL< I) 
8Ufn = <UI ( 1 )+Ul(2) ):*(X(2)-X( 1 ) >/2. 

SIJM2=(U.r <2)+WI(3) )>K(X<3)-X<2) )/2. 

5UM3=(UT<3)-1“UI(4) ))K(X(4)-X(3) )/2. 

SUM4= ( Ui T ( 4 ) +U1 1 ( 5 ) ) * ( X < 5 ) -X ( 4 ) ) /2 . 

USUM=SUm +5UM2+SUM3+SUM4 
IF(M-l) 428,426,428 
URF.Q=UCHAN/( AS)«(2S)<JRS) 

DIFF=Ai:(3(J4REQ-WSUrt) 

GO TO 430 

WR£Q=uJCHAH/( AR*7R*RR ) 

DIFF=ABS(UREQ-USUM) 

TAL=TOLi.YU'REQ 

IF (DIFF-TAL) 432,432,434 

VA(3)=VA<3) 

CODE=20 . 

GO TO 442 

IF(USUM-UREQ) 436,432,438 
CONTINUE 

IF(PRAT(i) .LT.PRATCR . AND.PRAT(5) .LT .PRATCR) GO TO 470 
^A(3)=VA(3)*(1 . 00*^OIFF/WREQ*i .01) 

GO TO 44? 



UA(3)=VA(3)*(i.00-DIFF/UREQ*l .01) 

WPER(1)=0. 

UipFR(2)=SUMi/l4SUM 

UP ER ( 3 ) = < SUM 1 +SUM2 ) / WSUM 

UPFR(4)=(SUMl+SUh2+SUM3)/USUM 

UPER (5)=1 . 0 

TF(IND^l) 450,423,423 

URtT£(6,422) (WI (I) ,1=1,5) 

FORMAT(/20H FLOW INTEGRAND 1^5, F10.5) 

UR I TE ( 6 , 424 ) SUM 1 , 5UM2 ; SUM3 , SUM4 , WSUrt 
FORMAT (/15H SUMS i-4, U3UM, 5Fl0.S) 

WRITE (6,440) WSUM , UREO , UA ( 3 ) 

FORMAT (35H REF FLOWS , COMPUTED-REQUIRED , AX UAL,?.F10.4, 
WRITE (6,444) WCHAN,WLBM 

FORMAT (/3<:»H REF FLOW RATE CHANNEL-SQUARE INCHES, F8. 5, 
-L.8M/SEC,F8,5) 

WRITE (6,446) M 

FORMAT (/30H STREAMLINE FLOW FRACTIONS, M=I2) 

WRITE (6,448) X(2) ,WPER(2) ,X(3) ,WPER(3) ,X(4) ,WPER(4) 
FORMAT (6F10 .4) 

GO TO 450 



FiO .2) 

18H FLOW RATE 



IN = 1 

RETURN 

END 



SUBROUTINE SLINE ( RR , X ,DUDX , WPER2 ,UPER 1 , HF ,U , DHEDX , S , DSDXl , 
)UARF,RRF ,FC1 ,FC2^,C0DE,M,B) 

DIMENSION RR( 10) ,X( 10) ,DU1)X( 1.0) ,UI ( 10 ) , WPER2( 10 ) , WPERl ( 10 ) ,HE( 10 ) , 
:(tDHEDX(i0) ,S( 10) ,DSDX1 (10) ,U< 1 0 ),B(20) 

COMMON/TOL/TOLl ,T0L2,T0L3,T0L4 

COMMON/IWl/IND, INZ, IWR 

N7=0 

SAUE=RR<3) 

CODE=l . 

DO 700 1=1,4 
J=I + i 

700 DWDX(I)=(WPER2( J)-WPER2(I) )/(X< J)-X(I) ) 

N=0 

DO 720 1=2,4 
K = I + 1 
J=I-1 

IF (ABS(WPER2(l)-UIPERi(I) )-T0L2) 716,716,702 
702 IF (WP£R2(X)-WPER1.(T) ) 704,716,708 
704 XN=X(3 ) + (WPER1(1)-UIPER2 (I))/Dw6x(J) 

IF(M-i) 706,712,706 
706 St.= (HE(K)-HE(l ) )/(X(K)-X(I) ) 

DEL=2.)i(SI.-DHEDX(I) )/(X(K )-X(I) ) 

DHEDX ( I ) = DHE 1)X ( 1 ) -»‘F)F1 % ( XN-X ( I ) > 

HE( I ;=HE( T )-»*nHEDX( T ))|c(XN-X( I ) ) 

S1. = (S(K )-S( I ) )/(X(K)-X(I) ) 

DEL=2. J<c(SI -DSDXl ( I ) ) / ( X ( K ) -X ( I ) ) 

DSDXl ( I ) =DSDX 1(1) -hDFF. * ( XN-X ( I ) ) 

S ( I ) =S ( I ) +DSDX1 ( I ) * ( XN-X ( I ) ) 

GO TO 71? 

XN=X( 1. )-(UPFR2( 1 )-WPFRl (I ) ) /DWnXd) 
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IF (h-1) 710,712,710 

SI = ( HE ( :r ) -HF (J))/(X(X)-X(J)) 

0EL=2.*(DHFPX( I)-SI.;/(X(I )-X( J) ) 

DHRI)X( 1 )=DHFDX( 1 ) +PFI )^ ( XN-X ( I ) ) 

HF( I )=HF( T)4-DH£DX( I ) ) 

SL=(5( 1 )-S< J) )/(XCl )-X( J) ) 

DEL =2 . * ( DSnXi ( I ) -SL) / ( X ( I ) -X ( J ) ) 
nSDXK I )=nSDXi(l)4.nFL4c(SN-X(I) ) 

S ( T ) =S ( X ) -♦•DSDXl ( I ) ( XN-X ( I ) ) 

RR(l) = XN)tcSAVF 
GO TO 718 
N=N+1 

IF (N-3) 720,730.720 
U(I)=U( J)*XN/X(1) 

CONTINUE 

DO 722 1=1,5 

X(I)=RR(I)/RR(3> 

Fni=RR(3)/SAUE 

FC2=FCi:*:ic2 

RRF=RR(3) 

ARF=B( 1 )+P(2)*RRF+B<3)*RRF*)|c2+B( 4):lcRRF*tv3-t-B(S)5KRRF**4 
IF(IND-l) 732>721^721 



!,^29 



XNEW HENEW 



OHEDX 



S-NEU DSDXi) 



IF (M-1) 729,73 
WRITE (6,724) 

FORMAT (/47H SLINE 
DO 728 T=i,5 

FORMAT ( I4'F9.4,F9.2,F9.4,F9. 6.F9.5) 

WRITE (6,726) I , X ( I ) , HE ( 1 5 , DHEDX ( I ) , 5 ( I ) , DSDXi < I ) 

GO TO 732 
GODE=40 . 

RETURN 
END 

SUBROUTINE ALOSl ( 7ETAS .ZETAPS) 
niMENSION ZETAS(IO) .7ETAPS(i0) 

COMMON/ Al 1 /ALFA! (10) ,Ui(10) , TTO , RPM , RS ( 1 0 ) ,SI ,TNX ,H ,D ,CI ,Ti ( 10 ) , 
*Pl(i0) ,TQ,TF1 jALI ,BF5P XX,ANG20,AMSl(in),S,TN,C,TE,AL.S0;THO, 
5|cCO,TFD,U(i 0 ) ,D11,D10,D2I ,D20,ANG2IiALFAX,TT ,1 .PTO , i^^LO 1 AMC 
COMMON/ A1.2/BETA2 U 0 ) , BETAl ( 1 0 ) , BFTAn( 1 0 > : 02 ( 1 6 ) , TTE ( 1 0 ) , U2 (10), 
*STR ,TNXR ,HR ,DZ,CIR ,TlPC,SZ,TNR,CR,BOR ,TN0R , COR , ALX R , ALR , ALOR , 
tP?. ( I 0 ) , WU2 ( 1 0 ) , W1 ( ( 0 ) , TEIR , TFR , TEOR , Di IR , D t OR , BETAZ , BETAI , ANM , 
XtTIR,lR;TnR,STAlll()0) 

CnMMUN/GAS/CP,GAM,EhME,ERRE,EXPl ,EXP2,VIS2,UIS3 
COMMON/CSS/CjjG.Qi 
nnMMON/IUX/IND,lNZ,IWR 
C0MM0N/AU5/XCL 
COMMON/ARE/REE 

COMMON/CnZ/ICOR ,IC07,I1NC,IAI,ICL,IAN,IC0N 
C0MM0N/ABA/BA17,BLEX , 

C0MM0N/TR5/TRA5 V * 

DO 6001 MACC=1,5,2 
TRA1=90 .-ALFA1(MACC)*57.29578 
TRA2=U1 (MACC)*.3048 
TRA3=TTO/l.a 

TRA4=RPM¥3, 14159/30 . :HRS ( MACC) /12 . » , 3048 
IF(MACC-3) 6002,6003,6004 

CALL TRAU? (TRAj,90. , TR A? , TR A3 , EMME ,GAM , SI , TNI , H , D , TR A5 , Cl , TR A4 , 
*0 . ,TRlto,TRA7,TRA6,TRA9,ZETAS(MACC) ) 

GO'tO 6001 , 

6003 CALL TRAU2 (TRAl ,90 . ,TRA2,TRA3 ,F.MME,GAM,S ,TN, H,D,TRA5,C ,TRA4, 
*0 . ,TR26,TRA7,TRA8,TRA9,Z£tA5(MACC) ) 

GO TO 6()01 

6004 CALL TRAU2 ( TR A 1 , 90 . , TR A2 , TR A3 , EMME , GAM , SO , TNO , H , D , TR A5 , CO , TR A4 , 
*0. ,TR36,TRA7,TRAfci,TRA9,Z£tAS(MACC) ) 

6001 CONTINUE 

TFdCOZ.LT.S) GO TO 2026 
DO 2027 MACC=t,5,2 
XY=ZETAS(MACC)/(i ,-7ETAS(MACC) ) 

ZETAS(MACC) = ( ( (1 .+XY)/(l.+XY*Pl(MACC)/PTn) ) *:1CEXP?-1 , )/ 

»( (PTO/Pl fMACC) )**FXP2-1 . ) 

IF ( IC07 , FO . 8 ) ZET AS ( MACC) = ( ( ( 1 . -hXY ) / ( 1 . +XY*BESP ) ) )K*EXP2-i . ) / 
*(1./BESP»)1^EXP2-1. ) 

2027 CONTINUE 
2026 CONTINUE 

XF(IC0N.NF.3) GO TO 31 
DO 32 1=1, 5, 2 
ZFTAPS( T)=ZFTAS(I) 

IF(XC0N.NF.2) GO TO 33 
TFdCOR .FQ.4) 7FTAPS ( 1 )=TR 16 
TFdCOR . F0.4) ZETAPS(3)=TR?6 
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IFdCOR .EQ.4) ZETAPS(5)=TR36 

33 CONTINUE 

rF(lCON.NE.l) GO TO 34 
r-n 35 1=1, 5. 2 
35 7FTAPS(I) = .53|cZETAS<I) 

34 CONTINUE 

DO 66 K=l,5 
66 CONTINUE 
RETURN 
END 

SUBROUTINE AL0S2 < ZFT AR , ZETAPR ) 

DT MANSION ZETAPSdO) ,ZETA3(1 0) .ZETAPR (10) .ZETAR(IO) 
COMMON/ALI/ALPAKIO ) ,U1(10) dTO , RPM , RS ( 1 0 > ,SI , TNT , H D ,CI , T 1 ( 1 0 ) , 
-1cPi(lO),Tn,TEI ,hLT ,BF^>P,XX,ANC20jAhS{ (10) ,s'tN,C.TF,Al.S01tN 0, 
*CO,TEO,U(iO) ,D1I ,Din . D21 . D20 , ANC2I > ALFAX > TT ,T ,Ptn, ALOI aMC 
CnrihON/AL2/BeTA2< 10) ;BETAi ( l6 ) .BETAfK 10 > , W2( iO ) , TtF ( 1 0 ) , U2 ( 1 0 ) , 
)KSIR ,TNXR .HR ,DZ,CIR .TXPC,SZ,TNR ,CR .50R ,TN0R .COR .ALJR >AI R ,ALOR . 
)KP2(10) ,1402(10) ,Wl 0.0), TtIR,TFR; T FOR, DiXR, DIOR, B£TA7,BETAI,ANA, 
)CTIR,TR,TOR,STALl(in) 

COhMON/UARi/RCdO) . R Sni.D2 . R S0LD3 > RSDI. D4 . aSFO . RSFO , PR FO . ARFO , 
*RR(10) .RR0LD2.RR0LD.3,RR0LD4,CU,CK,UA1 (10) , DALF ( 1 0 > , DBET ( 1 0 ) , 



)1cASF>AMS,Bi<205 
CnHMfjN/VAR2/I 



....... v^,,^../B6(20) ,7R .ZS.ARF ,B2(20) ,PR ,AMR 

CnNhON/UAR3/PTE(10) , R SI , RS3 , RS5 , T2 ( 1 0 ) 

C0MM0rVUAR4/BRl ,BR? >BR3,RR1 ,RR3,KR5,UA2(1Q) 

C0MM0N/ABA/BA17.BLEX 

CnrtrtON/GAS/nP,GAM,F.HME,ERRE,EXPl,EXP2,UIS2,VIS3 

CnhM0N/CSS7CJ,G,Qi 

COrtMON/IWI/IND,INZ,IWR 

COMMUN/AUR/XCL 

C0H?10N/ARE/REE 

COrtMON/COZ/TCOR , ICOZ , I INC , I AX , ICL, I AN , ICON 
C0hr10N/TR5/TRAS ' 

IFCICOR .NF,4) GO TO 60l0 
DO 6011 HACC=1,5,2 
TRA1=90 dBETA2(hACC))<c57. 29578 
TRAX=90 . -BFTA1(MACC):*57. 29578 
GXUD=BETAO ( MACC ) -BFT A1 (MACC ) 

IF ( 1 INC . EO . 1 . AND . GIUD . GE , 0 . ) TR AX=90 . -BETAO ( MACC ) %S7 . 29578 

TR A2=W2( hACC )* , 3048 

TRA3=TTE<HACC)/1 .8 

TRA4=U2 ( MACC )*, 3048 

XF(hACC-3) 6012,6013,6014 

6012 CALL TRAU2 (TRAl ,TRAX .TRA2 .TRA3 .EMME ,GAM .SIR ,X NIR ,HR ,DZ,TRA5,CIR , 
)*cTRA4.TIPC,TR16,TRA7,TRA8,TRA9,ZETAR(MACC) ) 

GO TO 6 0 1 i 

6013 CALL TRALI2 ( TR A i . TP AX , TP A2 . TR A3 , EHME , GAM , SZ , TNR , HR , DZ ,TR A5 , CR , 
:*:TRA4 .TIPC,TR26,TRA7,TRA8,YCL , ZETAR < MACC ) ) 

GO X6 60li 

6014 CALL TRAU2 (TRAl . TRAX ,TRA2 .TPA3 .FMME .GAM ,SOR ,TNOR ,HR ,DZ ,TRAS,COR , 
*TRA4.TIPC, TR36,TRA7,TI^A8,Tf<A9, ZETAR (MACC) ) 

6011 CONTINUE 

DHl=CP5»cTT05<<( 1.-(P2(3)/PTO)^^<5 |cRXP2) 

PS T = 1 . / ( i . -YCLJKDH 1*G*C J/U2 ( 3 ) /WU2 ( 3 ) ) 

ZEZE=ZETAR<3) 

DO 6015 MACC=1,5.2 

ZETAR ( MACC) =ZETAP < MACC) + ( 1 .-ZEZE)ilc(l .-PSI^^cPSI) 

6015 CONTINUE 
6010 CONTINUE 

IFdCOZ.LT.S) GO TO 2046 
DO 2U47 MACC=1,5 
XY=ZETAR(MACC)/(1 . -ZETAR < MACC ) ) 

ZETAK(MACC) = ( < (1 .-hXY)/( 1 . +XY*P ? < MACC ) /PTE ( MACC ) > )*)KEXP2-1 . )/ 
(PTE(MACC)/P2(MACC) ))*C)*(EXP2-1 . ) 

IF(ICOZ.EO.R) ZETAR < MACC ) = < ( (1 .+XY)/(1 .+XY*BESP) ))*t*EXP2-l . )/ 

*( 1 . /BESP**EXP2-i . > 



>046 CONTINUE 
?047 CONTINUE 

IFdCON ,NE,3) 
30 DO 32 1=1,5 2 



GO TO 31 

32 7ETAPR( J )=ZI?:TARd) 

31 IFdCON.NE.2) GO TO 

TFdCOR .FQ.4) 
IFdC0R.Ef;!.4) 

XFdCOR .EQ.4) 

33 CONTINUE 
IFdCON.NE. 1) GO TO 34 
DO 35 X=i,5,2 

35 ZF.TAPRd)=.S^ZLTAR(I) 



33 

7.ETAPR<1)=TR16 
ZETAPR(3)=TR26 
ZE1 APR(5)=TR36 
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34 CONTINUE 

. 7F:TAR(2)= ZFTARCl ) + ( RR ( 2> -PR 1 ) / ( RR3-RR i > ^ ( 7E TAR ( 3 ) -7ETAR ( 1) ) 
ZETAR(4) =7ETAR ( 3 ) + ( RR ( 4 ) -RR3 > / ^ R R 5-RR3 ) ( ZETAR ( 5 > -ZETAR ( 3 ) ) 

7ETAPR(2)=7ETAPR( 1 ) + (RR(2)-RR1 )/<RR3-RRl ) 5K < ZET APR ( 3 ) -ZF.l APR ( 1 ) ) 
ZETAPR (4)=7ETAPR (3)+(RR (4)-RR3)/(RR5-RR3)JK(ZETAPR(5)-ZETAPR (3) ) 
DO 7001 1=1,5 
7001 CONTINUE 
RETURN 
END 
C 

FUNCTION MAURA(TH,TE,SP) 

C TH=THROAT OPENING 

C TF=TRAILJNG EDGE THICKNESS 

C SP=BLADE SPACING 

ARG1=TH/SP 

TFRMl=ATAN(SORT(i . - ( ARGi>*ciK2 ) )/ARGi) 

TFRn2=TFRr11«l80 ./3. 14159 
TFRM3=i . -TFRh2/90 . 

TF.RM4= ( 4 . .yTE/SP ) *TERM3 
ARG2=n H/SP ) +TFRH4 

UAURA=ATAN(SQRT(1 ,-(ARG2**2) )/ARG2) 

RETURN 
END 



C 

C 

C 



1 



FUNCTION XPO( ANG£ ,ANG2) 

C0hh0N/TRA/XP0i(5,8) ,XP02(6,B) , ALFl ( 8 ) , ALFOi (5) ,ALF02(6) , 

:«cY( 10) . Y1 (10) ,Q(6) ,RX(30 ),RY(30),IR(30),Z(6),Ci(4,8),C2(4,8) 
TF(ANG2-80.) 1>2,3 
CONTINUE 
DO 4 1=1,8 
DO 5 J=1.4 

5 Q(J)=C1(J.I) 

4 Y(I)=YC(ANG2,Q,3) 

GO TO 10 

2 CONTINUE 
DO 6 1=1,8 

6 Y(I)=XP0i(5,I> 

GO TO 10 

3 CONTINUE 
DO 7 T=i,8 
DO 8 J=1^3 

B 0(J)=C2(J»I) 

7 Y(I)=YC(ANG2,Q,2) 

10 CONTINUE 

DO 11 1=1,7 

IFCANGi.GE.ALFKI) ,AND.ANGl.LE.ALFi(I + i)) GOTO 100 
IF(ANGl.LT.ALFKl)) GO TO 101 
IFCANGl ,GT. ALFKS) > GO TO 102 

11 CONTINUE 

100 CONTINUE 

XPO=Y(I )+(Y(l+l)-Y(I> )/(ALFl( I+l )-ALFl(I) ) * ( ANGl-ALFl ( I ) ) 
IF(ANG2. LT. 40 ) XPO=0 , 09-(0 . 09-(XP01 ( 1 , I )+XPOi ( 1 , I+i ) )/2. )* 
3 |c(ANG2-20 . )/20. 

RETURN 

101 XPO=Y(l) 

RETURN 

102 XP0=Y(8) 

RETURN 

END 

FUNCTION nSIM(Ul ,TO,FhhE,GAM) 

FRRE=848 . *9 , S0665/FHME 

AST=SQR T ( 2 , >KGAM / ( G AM+ 1 . ) ^KER R E*TO ) 

AMACH=U1/AST 
■ “ .8) CSIM=i. 

,8) GO TO 1000 

1) CSlh=i .-0 .22/0 ,3)lc( AhACH-0 .8) 

IFCArtACH.LT . 1 . 2 . AND , AhACH . GT , 1 , t ) CSIh=0 .78+0 . 15/0 . 1 * ^MACH-1 ,i: 
IF(AMACH.GT. 1.2) CSIN=0 . 92+1 . 5/ . 2* ( ArtACH-1 .2) 

1000 RETURN 
END 

SUBROUTINE CID ( ANGl , T. DFL ,CSID , PSID ,PSXF , Hh , DM) 

DIMENSION X(7),Yi(7),y2(7> 

FF=1 .-DFIVT/SIN(ANG1) 

DATA X/ 15. ,20 . ,25. ,30 . ,45. ,60 . ,90./ 

DATA Yl/ i! 06,1 .1 .1 .:l7,1.225,i:63,2.i ,2.45/ 

DATA Y2/n 01 6 . 0 . 031 5 , 0 . 049 . n 0-70 ^ ^ , 1^6 0.260,0,4 / 



IFCAMACH.LE, 0 . 
IF(AMACH.L.E. 0 . 
IF(AMACH.LF. 1, 
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C 

C 

C 

C 

C 

C 

c 

c 

n 

c 

c 

c 

c 

c 

c 

c 



Y=+(Yi(I+i)-Yi<I) )/<X(I+i)-X<I) )# 



(Y2< I+i)-Y2(l) )/<X(I+l)-X( I) )* 



A=.ANGi)<cl80 . /3. 1415 
. i:-0 i 1 = 1 » 6 

IF(A.LE .X( 1 ) ) Y=i .+0 . 06*A/1S. 

IFCA.GE.XCI) .AND. A.LE .X(I+i) ) 

)K(A-X(I) )+Yl(I) 
lECA.L e ,X( t ) ) Z=Y2<1)*A/X<1) 

IFCA.GE.Xd). AND . A . LE . X < I ->*1 ) ) 

*<A-X(I))+Y2(I) 

L CONTINUE 

IF(A.GT.X<7> ) Y=l. 

IF(A.GT.X<7) > Z=t . 

CSID=i . -i-<Y-l . )*2.*(1 .-EF) 

PSll>-Z%A .tiX ,-F>))K<l .-EF) 

P5.TF=0 . 025/0 . 09)#: Hh*HM/DM/DH 

RETURN 

END 

FUNCTION CSIUCXPO.CBIP .T,ANG1,AH) 

C3IU=XP0*CSIP*T*SIN(ANG1)/AH 

RETURN 

END 

FUNCTION CSIRCS.AH.Ul ,ANGi,UM,XP) 

SL=S/AH 

I F ( SL . LE . 0 . 4 ) XL=XP a . 65/ . 4*SL 

IF(SL,GT,0 .4.AND.SL.LF.0 .8) XL=XP*(0 .65+0 .45/0 . 4ik ( SL- . 4 ) ) 
TF< SL . GT . 0 . 8 . AND . SL . I E . t . 5 ) XL=XPt< 1 . 1 + 0 . 04/0 . 7)t{ (SL-0 . 8 ) ) 
IFCSL.GT. 1 .5) XL-XP>K( 1 .5+0 .6/1 ,71t(SL-l .5) ) 

ASC=Ui)KSTN<ANGl )/UM 

XPO=0 . 025+0 . Oib/0 ,64:<(ASC*ASC 

JFCASC.LT. .2) XRO=.024 

IF<ASC.GT. .95) XRO=.0475 

CSIR=XRO»XL 

RETURN 

END 

FUNCTION ALEAK < DELPFT, DM >AL,CLE, ALFA!) 

Ci = 0 . 82-0 . 075JKDELPET 

ALEAK=Ci¥<DM+AL):*cCLE/DM/AL/COS<ALFAl) 

RETURN 

END 



SUBROUTINE CHBFT ( X . Y . N . A . M , RX , RH , R ) 

DESCRIPTION OF PARAnFTERS: 

X ARRAY OF ABSCISSAE DIMENSIONED RFAL^A X(N) 

Y ARRAY OF ORDINATES DIMENSIONED RFAL*4 Y<N) 

N NUMBER OF SAMPLE POINTS (INTEGER) 

A ARRAY OF I HE OUTPUTTED POLYNOMIAL COEFFICIENTS 
DIMENSIONED AT L.EAST ACM+2) (REAL*<)) 

M ORDER OF DESIRED APPROXIMATING POLYNOMIAL 
RX WORK ARRAY DIMENSIONED AT LEAST REAI.*4 RXCM+2) 

* RH UORK ARRAY DIMENSIONED AT LEAS! RfALtA RH<M+2.) 

R INTEGER UORK ARRAY DINENSIONED AT LEAST R<M+2) 

-NOTE: DIUIDED DIFFERENCES AND NEUTON^S INTERPOLATING FORMULA IS 
USED FOR COMPUTING THE POLYNOMIAL COEFFICIENTS . 

REAL NEXTHI 
INTEGER RI,RJ.R<i) 

DIMENSION X<i)>Y<i),A(l),RX(i),RH<l) 

MPLUS1=M+1 
hPLUS2=M+2 
PREUH=n . 0 

DETERMINE INDEX UECTOR FOR INITIAL REFERENCE SET 

R(i)=l 

R(MPLUS2)=N 

D=(N-i)/MPLUSi 

H=D 

DO 1 I=2,MPLUSi 
R(I)=H+1.0 

1 H=H+D 

2 H=-l . 0 

SELECT M+2 REFERENCE PAIRS AND SET ALTERNATIUE DEVIATION VECTOR 

DO 3 I=i,MPLUS2 

RJ=R<I) 

RX<I)=X<RI) 

A(I)=Y<RI) 

3 r"S7V,=h , 
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C 

C 



C 

C 



C 

c 

c 



12 
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11 



110 



14 

15 



COMPUTE M-M LEADING DIVIDED DIFFERENCES 
DO 4 J=i,MPLUSl 
I1=MPLUS2 
AI1=A( II ) 

RHT1=RH( ID 
T=nPLlJ81 

DENOM=RX( II )-RX(I-J+l) 

AI=A(I) 

RHT=RH(T ) 

A( U )=(AT1-AI )/DENOM 
RH(11)=(RHX1-RHI )/DENOM 
11 = 1 
AI1=AI 
RHU=RHI 
1 = 1-1 

TF(I-J) 4,5,5 
CONTINUE 

EQUATE (M-M) THE DIFFERENCE TO ZERO 10 DETERMINE H 
H=-A ( MPLUS2 ) /R H ( MPLU82 ) 

WITH H KNOWN, COMBINE THE FUNCTION AND DEVIATION DIFFERENCES 

DO 6 X = 1. .MPLUS2 

A(I)=A(I^+RH(I)HH 

COMPUTE POLYNOMIAL COEFFICIENTS 

J=M 

XJ=RX(J) 

I = J 

AI=A(I) 

JPLUS1=J+1 

DO 3 11=JPLUS1,MPLUS1 

AI1=A(ID 

A<I)=AI-XJ*AI1 

AI=AI1 

T = I1 

J=J-i 

TF( J-1 ) 9,7,7 
CONTINUE 

IF THE REFERENCE DEVIATION IS NOT INCREASING MONOTONICALLY 

THEN EXIT 

HMAX=ABS(H) 

IF(HMAX .GT.PREVH) GO TO 29 

A(MPLUS2) =-HMAX 

RETURN 

FIND THE INDEX. IMAX, AND VALUE, HIMAX, OF THE LARGEST ABSOLUTE 

ERROR FOR ALL SAMPLE POINTS 

A<MPLUS2)=HMAX 

PREVH=HMAX 

IMAX=R(1) 

HIMAXi^H 

J=1 

R,T=R<J) 

DO 110 1=1, N 
IF(T.EQ.RJ) GO TO 11 
XI=X(I) 

HI=A(MPLUS1) 

K = M 

HI=HI*XI+A(K) " 

K =K *“ 1 

IF(K-l) 112,12,12 
HI=Y(I) 

ABSHI=ABS(HI) 

IF(ABSHI.LE.HMAX) GO TO 11 

HMAX=ABSKI 

HIMAX=HI ' 

TMAX=I 

GO TO 110 • 

IF( J.GE.MPLUS2) GO TO 110 

J=J+1 

R.T=RCJ) 

CONTINUE 

IF IHE MAXIMUM ERROR OCCURS AT A NGNRFFFRENCE 
POINT UT7H THE NEARFSF REFERENCE POINT HAVING 
SAME SIGN AND REPEAT 
TFdMAX.EQ.Rd) ) RETURN 
DO 14 1==2.MPLUS2 - ' 

TF(IMAX.LT.Ra)) GO TO 15 
CONTINUE 
I=MPLUS2 : 

NEXTHI=H 

TF((l-T/2*2) .NE.O) HFXTHI=-H 
IXHXMAX^NFXTHI.GF.H) GO TO 115 



POINT^ EXCHANGE THIS 
AN ERROR OF THg 
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0959 

0960 

096i 

09<S2 

0963 

0964 

0965 

0966 

0967 

0968 

0969 

0970 

0971 

0972 

0973 

0974 

0975 

0976 

0977 

0978 

0979 

0980 

0981 

0982 

0983 

0984 

0985 

0986 

0987 

0988 

0989 

0990 

0991 

0992 

0993 

0994 

0995 

0996 

0997 

0998 

0 999 

1000 

1001 

1002 

1003 

1004 

1005 

1006 

1007 

10 08 

1009 

1010 

1011 

1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 

1021 

1022 

1023 

1024 

1025 

1026 

1027 

1028 

1029 



C 

C 



IF(IMAX.r,E.R(l)) GO TO 116 
Jl=nPLUS2 
' J=f1 . 

117 R(J1)=R(J) 

J1=J 

J=J-1 

TF(J-l) 118,117,117 

118 R<1)=IMAX 
GO TO 2 

116 IF(lMAX.LF.R(hPLUS2)) GO TO 120 
J=1 

no 121 J1=1,MPLUS2 
R(J)=R(J1) 

121 J=J1 

R(MPLUS2)=IHAX 
GO TO 2 
115 R(T)=IrtAX 
GO TO 2 

120 R<I-i)=IHAX 
GO TO 2 
END 

SUBROUTINE TRAUl 

COrthON/TRA/XPQl (5,8) • XP02(6,B) , ALFl (8 ) ,ALF0i(5) ,ALF02(6) , 

*Y(i0) ,YK10) ,Q(6) ,RX(30),RYt30) ,IR(30 ) ,2(6) ,C1 (4,8) ,C2(4,8) 

DO 6 1=1,8 
DO 7 J=l,5 

7 Y(J )=XP01 ( J,I) 

DO 8 J=l,6 

8 Y1(J)=XP02( J,l) 

CALL CHBFT ( ALFO 1 ,Y,5,Q,3.RX,RY,IR) 

CALL CHBFK ALP 02, Yi ,6,2,3, RX,RY,IR) 

DO 12 J=l,4 
C1(J,I)=Q( J) 

12 C2(J,I)=Z(J) 

6 CONTINUE 
RETURN 
END • 

SUBROUTINE TRAU2 ( ANGi , ANGO ,U1 ,TO ,ErtME ,GAh , T ,DEZ ,HM , DM ,CSIP ,S ,Urt , 
)*<CL>RPR0,R2,R3,YCL,RT0T) 

CSip=i. . 

R=XPO(ANG1. ,ANGO) 



P1=CSIM(U1 >TO,PMME,GAM) 
ANGZ=AN'G1*3. 1415/180 . 



10 



CALL CID( ANG7 .T ,I>F7 ,CSTD,PSID,PSIF,HM,DM) 
R2==CS lU ( R , CSIP , t , ANGZ , HM ) 

R3=CSIR (S,HM,Vi,ANQ7,UM,CSIP) 
RPRO=R*CSTP)KRl*CSID+PSIF-»'PSID 
IF(CL.LE.0. ) YCL=0. 

IF(CL.LE. 0 . ) GO TO 1000 

DEL=3. 1416-( ANGO+ANGl )*3, 1416/180 . 

ALF1=1 .57Q8-ANG2 
YCL=ALF AK ( DEL .DM , HM , CL , ALF 1 ) 
on RT0T=RPR0+R2+R3 

RETURN 

END 

FUNCTION YC(XBAR,Q,M) 

DIMENSION Q(6) 

YC=0 

IF(XBAR.EQ,0.) YC=g(i) 

IF(XBAR ,EQ. 0. ) GO TO 10 



M1=M+1 

DO 1 1=1, Ml 

1 YC=YC+Q(I))*tXBAR)|c«(I-l) 
10 CONTINUE 
1000 RETURN 
END 
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.^.SHORT T = 00003 IS ON CR0002S USING 00024 BLKS R=0000 



0001 

0002 

0003 

0004 

0005 

ooos 

0007 

0008 
00 09 
0010 
0011 
0012 
00 5 3 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

mt 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 
0 063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 



FTN4>L 



C 

C 



C 

C 

C 

C 

c 

c 

c 



PROGRAM BHORT(S) 

DIMENSION INAM(3) 

DIMENSION NAME(3) 

CnMM0N/ABA/BA17 ,BLEX 
COhhON/CUR/COSL (10) 

COMMON/TOL/TOLl , T0L2 , T0L3 , T0L4 
C0MM0N/TRS/TRA5 

r.OMMCiN/f'AB/CP,GAM.F.MME.ERRF-,F.XPi ,EXP2,UIS2,VIS3 

C0Mn0N/C02/IC0R , IC07, IINC, lAl , ICIl , IAN , ICON 

COhMON/hAC/IN 

COMMON/IUT/IND,INZ, IWR 

r.OHMON/AUS/XCL 

COrtriON/CS8/CJ.G»Ql 

COhhON/UARl/RC( 10 ) ,PB0LD2.RSnLD3>RB0LD4, ABFO.RSFO.RRFO, ARFO, 

*RR ( 10 ) ,RROLD2,RROLD3,RROLD4 .CV,CK ,UA1 ( lO) ,DALF( lo5 ,DBET( 10 ) > 
:1JABF,AMS,P1(20) 

COMMON/UAR2/B6(20 ) s ZR . ZS^ARF , B2 ( 20 > , PR > AMR , UUl ( 1 0 ) 
COMMON/UAR3/PTE( 105 , RBl , RS3 , R B5 , T2 ( i 0 ) 

C0MM0N/UAR4/BR1 ,BR?,BR3>RR5 , RR3 , RR5 , UA2 < 1 0 > 

mnnON/UARS/PRATI (lA) ,R1NCI<10> .ALFAIKIO) ,BETA11(10) ,ZETAl(i0), 
)K()2 ( 1 0 ) . Ai .FA22 ( 1 0 , BET ^22 ( 1 0 ) 

COMMON/UAR6/PT2( 1 0 ) ,TT2(10 ) ,PT1(10) ,DEUH(10) ,ALFA2( 10 ) ,UU2(10) , 
)#CWR2(10) ,T2S(10) ,T2IS(10) 

C0MM0N/UAR7/TTIS(5) ,BF.TAT(5) , ETAT (5 ) , ETAI ( 1 0 ) ,ETAS( 10 ) ,ETAR( 10) , 
*RSTAR ( 1 0 ) . AK I S ( 1 0 ) , PS TR ( 1 0 ) 

C0MM0N/UAR8/DR1 (10 ) .AMUl(lO) ,AM'J2(10) .BFTFTd 0) .PRATlTdO) ,AMR2( 
:*:10> , YS( 10 ) ,X1 (10 ) ,AfiFAl (10) ,7FTAPS(10 ) ,UPER1 (10 ) » YR ( 1 0 ) , X2 ( 1 0 ) 
COMMON/UAR9/ZETAR (5 0 ) .ZETAPR (10) , AB ( 1 0 ) , AR ( 1 0 ) , S 1 1 (1 0 ) ,SI2d0) , 
)KSldO) .DBDXldO) ^UJTKIO) ,HEdO) 

CnMMGN/UARlO/WUl (5 0) >DHEDX(5 0) ,DBDX2( 5 0) ,RIidO),RI2dO), 

)itR 13 (1 0 ) , R 14 ( 1 0 > , R T ( 1 0 ) , SR 1 d 0 ) , BR 2 ( 1 0 ) 
rOMMON/UARll/YOLDdO) ,AA(10) .BRdO) ,PRAT2dO) ,UPER2dO) > 
*DWDX(10) .TtlSdO) ,PRAT3(10) . BSd 0 ) , ALFA d 0 ) 

COMMON/U AR 12/bEI A ( 5 0 ) , DELR (1 0 ) , UPERO ( 1 0 ) , ZFTAS (10), ZETAR 1(20), 
>K7.ETAR3(20 ) , ZETAR5 ( 20 ) , R 1 ( 20 ) , A1 ( 20 ) , T 1 0 ( 20 ) 

CnMrtON/UARi3/Sri (20 ) , IRP ( 20 ) , R2 ( 20 ) . A2(20) ,R INC ( 20 ) ,DR ( iO ) , 
)KBFTO(10) .STALIKIO) .AREA2(iO) ,UR1(10) 
C0MM0N/UAR14/ULBM,PRAT5,0MFG 

COMMON/ALI/AlFAI (10) ,UldO) ,TTO,RPM,RS(10) ,SI,TNI sH >D,CI ,TldO) , 
)iJP I ( 1 0 ) , TO , TE I , AL I , BESP , XX , ANG2Q , AMS 1 d 0 ) , B , TN , C , TE , AL , SO , TNO , 

3»(SIR,1NIR .HR ,DZ,r:iR ,TIPC.SZ,TNR'CR .BOR ,TNOR .COR ,AJ IR ,ALR JaLOR . ' 

*P2(iO) ,wfl2(10 ) ,W1 (iO) ,T£IR,TER,TEOR,D{lR,DiOR,6ETAZ^BETAT,AN«, 
*TIR.TR,TnR,STALI(10) 

COMMON>ARE/REE 

COMMON/TRA/XPOi (5,8) .XP02(6,8) , ALFl (8 ) , ALFOl ( 5 ) ,ALF02(6), 

)KYdO) . Y1 (10) .Q(6) ,RX(30 ) ,RY(30), IR(30) ,Z(6) ,C1 (4,8) ,C2(4,8) 

DATA iNArt /28SH,2H0R,2HT / 

DATA NAME /2HP A , 2HRT , 2H2 / 

CALL TRAUl 
XX=1,25 
BLEX=^0.t5 . 

XCL=1 .35 

BESP = (i .+( GAM-1 . ) /2. 64) *5K( -GAM/ (GAM-1 . )) 



C 0 N D I T 1 0 N S * * * * * )|c)*c * )*c)tc)tc * 



***)|c**)lC}|c)H*OPERATING 
PT0=38.22 
TT0=626. 18 
RPM=12000 . 

PR=2.6 

C 

C 

c )|c***)<c)»:*)tclNITIAL APPROXIMATIONS5lc«*)|c)tc*)tc)|c)tc)iC)|o|c*)lc5«t 

AMC=.2247 
AMS*. 9 
AMR*. 7 

UAi(3)=26? .58 
^)A2(3)=262.58 
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0079 

OOSO 

0091 

0 082 

008 ^ 

0094 

0085 

008S 

0 087 

0088 

0 089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

U1S8 



4c ¥ 4c * * * 'ifc )(c * iK * .4C )(ci( )H * * * 4? * 

*4c4:5t4c4cSPfc.CI AL INPUT DAT A4c4c4c4c4c4c4c4c45 4c¥4c^i4c4c4c4c4c4c4c4c 
TOLl= . 01 
TnL2=. 01 
T0L3=. 01 
T0L4=. 01 

TND=0 
INZ = 0 
I14R=0 
IC0R=4 
IAI = 0 

TAN=2 ’ * 

ICL = 0 

TINC=i 

XCOZ=6 

IC0N=3 

4( 4C 4c 4^ 4( 4C 4( 4: 4C 4c 40K 4( 4c 4( 4( 4( « 4( 4( 4C 4c 4( 4c 4( 4C 4c 4( 4( 4c 4( 4C 4c 4C 4( 4( 4( 4c 4( 4c 4( 
4C4c4c4c4c4c4c4c4(TliRPTNE GEOMETRY 4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c4c 



A1 


(1) 


= 


.2126 


A1 


(2) 


= 


.22145 


A1 


(3) 


= 


.23035 


A1 


(4) 


= 


.23925 


A1 


(5) 


= 


.24815 


Al 


(6) 


s= 


.25705 


A1 


(7) 


= 


.26595 


Al 


(8) 


= 


.2.7485 


Al 


<9) 


= 


,28375 


Al 


(10 


) = 


=.2926 


A2 


(1) 


= 


. 1912 


A2 


(2) 


= 


.2.0305 


A2 


(3) 




.21495 


A2 


(4) 


= 


.22685 


A2 


(5) 


= 


.23875 


A2 


(6) 


= 


.25065 


A? 


(7) 


= 


.26255 


A2 


(8) 


= 


.27445 


A2 


(9) 


= 


.28635 


A2 


(10 


) 


=.2983 


AL 


= 1. 


088 


AUT = 1 


, 


088 


ALD=i 




088 


c= 


1.0 


03 


Cl 


= 1 . 


0 


03 


CO 


=1. 


0 


03 



C 

C 



C 

C 

c 



E=2.a065 

£1=2.8065 

EQ=2.806S^ 

T=.2252 

TI=,2252 

T0=.22S2 

TE=.03 
TEI=. 03 
TEO=.03 

TN=.0186 
TNI= .0186 
TNO=. 0186 

ALR=1 . 088 
ALIR=i . 088 
ALOR=l. 088 

CR=1 .003 
CIR=i. 003 
COR=l ,003 

FR=2.45 

£TR=2.45 



255 



0159 

0160 
0161 
0162 
0163 
016 ^ 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 
0201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 



E0R=2. 45 



TR=.2252 
TIR=.2252 
T0R=i.2252 
C 

te:r=. 03 
TEIR=. 03 
TEOR=. 03 
C 

TNR=. 0186 
TNIR=. 0186 
TNOR=. 0186 
C 

RC(1 )=2.764 
RC(5)=3.627 
C 

RS(1)=2,764 
RS(5)=3.627 
D 

RR(1)=2.693 
RR(5)=3.837 
D 

CV=.885 
CK=5, 0 
TIPC=.0i 
7S=31 . 

ZR=32. 

C 

RC(2)=SraRT(RC(5)lcRf:(5)/4.+3. /4.*RC(i)*RC(i) ) 
RC(3)=8nRT(RC<5))*cRC(5)/2. -»-RC< 1 ))KRC(i )/2 . ) 

RC(4)=5QRT(3 . /4 . :»cRn ( 5 ) *RC (5 ) +RG ( 1) *RC ( 1 ) /4 . ) 
RS(2)=Sf3RT(RS(5)*RS(5)/4.+RS<l )*RS(l)/4. *3. ) 
RS(3)=(RS(1)+RS(5) )/2. 

R5<4)=SQRT( < ( RS ( 1 ) **2 ) /4 . )) 

PR ( 2 ) =Sr;R T ( R R < 5 ) *RR ( 5 ) /4 , +3 . /4 . )i{RR ( 1 ) *RR ( 1 ) ) 
RR(3)=(RR(l)+RR(5) )/2. 

RR(4)=S0RT(RR (5))ICRR (5>¥3./4. +RR(1))*CRR(1)/4. ) 

DO 3300 I=l>10 
A=I 

R1(I)=RS(1 )+(A-i . )/9 
R2(I)=RR(1)+(A-1 . )/9 
3300 CONTINUE 
OI=Al(i) 

no=Al(iO) i 

0IR=A2(1) 

DOR=A2(10) \ 

DO 3711 1=1,10 

IF(RS(3) .LE.RKD) GO TO 3712 
CONTINUE 
CONTINUE 
T = T-1 

D=A1(I)4-(A1 (Ii-i)-Al ( I) )*(RS(3)- 
DO 3713 1=1,10 

TF(RR(3) .LE.R2(I>) GO TO 3714 
CONTINUE 
CONTINUE 
1 = 1-1 

nR=A2(I)+(A2(I+l)-A2(I) )*(RR(3)-R2(I) ) /( R2 ( I+i) -R2( I ) ) 
H=RS(5)-RS(1 ) 



,*<RS<5)-R8(1>) 

*(RR(5)-RR<1)) 



3711 

5712 



3713 

3714 



-R1(I))/(R1(I+1)-R1(I) ) 



C 

C 

C 



D=2.*RS(3) 

8=2. *3. 1416*RS(3)/2S 
SI=2.*3. 14165 *RS(1)/ZS 
50=2. *3. 1416)i<RS(5)/ZS 
57=2. *3. 1416*RK(3)/ZR 
SIR=2.^3. 1416*RS( 1 )/ZR 
S0R=2.*3. 1416*RS(5)/ZR 
DZ=2. *RR (3) 

HR=RR(5>-RR(i> 

STATOR OUTLET ANGLES BY UAURA rtETHOD 

ALFA! (3)=UAURA(0,TN.S) 

ANG2I =UAURA(OI ,TNI ,SI ) 

ANG20 =UAUR A ( DO , TNO , SO) 



DALF(i)=ANG2I- 
DALr (3:> = 0 . 



ALFA! (3) 
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0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 

0308 

0309 

0310 

0311 

0312 

0313 

0314 



C 

c 

c 



c 

c 



DALF(5)=ANG2Q-ALFA1 (3) 

ROTOR OUTLET ANGLE5 BY VAURA METHOD 

PETA2(3)=-1 . #OAORA(nR ,TNR ,S7) 

RETAI =-t .^nUAURACnTR.TNTR.SIR) 

BETAZ =-l . )fcOAOR A(OOR ,TNOR ,SOR ) 

DBETC 1)=BETAI-B£TA2(3) 

DBET(3)=0. 

DBET ( 5 ) =BET AZ-BETA2 ( 3 ) 

Xl(I)=RS( t)/RS(3) 

10 X2<I)=RR(I)/RR(3) 

RR0LD2=RR(2) 

RR0LD3=RR(3) 

RR0LD4=RR(4) 

RS0LD2=RS(2) 

R.S0LD3=RS(3) 

R5nLD4=RS<4> 

RS1=RS(1) 

R33=RS(3) 

RS5=RS<5) 

RR1=RR(1) 

RR3=RR<3) 

^ ( 5 ) 

CALL CHBFT(Rl,Af, j 0,B1 ,4,T10 ,8T1,1RR) 

CALL CHBFT(R2'A2Uoj B2!4,Tt0 ,5T1 , IRR ) 

ASF=Bi ( 1 )4-Bl (2)*RS(3)+B1(3):«RS(3))*C)K24-B1 (4>)KRS(3).ioK3 + Bl(5))KRS(3>*)K4 

ARF = B2( i)-fB2(2)*RR(3)+B2(3))»cRR(3)*)<c2+B2(4)5*{RR(3))*c*3 + B2(5)*RR(3)^)*c4 

ASF=.2S26 

ARF=.2447 

RSF-RS(3) 

ASFO^ASF 

RSFO=RSF 

RRF=RR(3) 

RRFO^RRF 

ARFO=ARF 

: INPUT PRINTING 



R1 



A1 



WRITE(6.67i) 

671 FORMAT(iHl,//60X,j2HINPUT PR INT3//40X >50H 
* ^R2 A2 /) 

DO 72 1=1,10 
WRITE <6, 73) Rl(I) ,A1(I) ,R2(I) .A2(I) 

73 FORMAT<40X,F10 .3,f^l0.4>i0X,Fl 0.3,F10.4> 

72 CONTINUE 

WRITE(6%74) ZS 7R TIPC CU CK 

74 FORMAT (/4SX,25HNUMBER "6f ^TATOR BL ADES = . F8 . 0/45X , 25HNUMBFR OF ROT 

.YOR BLADES = . F8 . 0/45X , 25HR0T0R TIP CI..FAPANCE = . F8 . 4/45X , 25HAXI 

^AL DISTANCE L = , F 8 . 2/45X , 25HCURUATURE FACTOR K =,F8.2/ 

*//55X, 16HBLADING GEOMETRY/) 

WR17E(6,75) 

75 FORMAK//30X,70H C E T TE TN 

^ A) R / ) 

WRITE(6,76)CI,EI,TI .TEI^TNI ,ALI,RS(1) , C , E , T >TE , TN , AL , RS ( 3 ) , CO , EO , 



*TO,TEO,tNn,AL6.R^<55 

76 FORMAT(30X,7Fl0.4/22Xv6H5TATnR,2X>7Fj 0 . 4/30X , 7F 1 0 . 4/ ) 

77 FORMAT ( 30 X : 7F1 0 . 4/22X ) 6HR0TnR :2X,7F10 , 4/30X, 7F10 . 4/) 

WRITE(6 ,77^C1R,E1R .TIP ,TE1R , TNt R , 4LIR , RR ( 1 ) ,Cr ,ER,TR ,TER .TNR ^ALR, 

*RR(3) .COR ,EOR ,TOR ,TE0R ,TNOR ,AL0R,RR(5) 

URITE(6,7ft) 

78 FORMAT(/40X,52HALL DIMENSIONS INDICATED IN THIS TABLE ARE IN INCHE 
*S/> 

WRITE(6,79)TC0R ,IAI .TAN.lCnZ, TINC.ICL,irON 

79 FORMAT(////,40X,27HCORRELATION SYSIFM ICOR = , I5/61X , 6HT AI =>15/ 
:»J6iX.6HTAN = , 1 5>6 1 X , 6HIC0Z =, 15/ , 61X , 6HT INC =, T5/61X ,6HICL =>I5/6 
*1X.6hICON =>I5) 

URllE(6,Bl )CP .EMME.GAM,UIS2,UIS3 

81 FORMAT(/20X, 91HGAS PROPF_RTIES CP MOL FClIl AR MASS 

GAMM OtSCOSTTY (1) UISCOSITY ( 2 ) /38X , 1 OH ( BTU/LB F);32X, 13 

*H(LBM /SEC FT) ,4X,13H(LBM /SEC FT)//36X,F9 . 3,SX,F1 0 . 3,9X,F7 . 3, 2E15 

^,Z/) 

CALL EXEC<8,NAME) 

END 
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H1N4 >L 

PROGRAM PART2(5) 

DIMENSION NAME(3) 

DIMENSION NAMR(3) 

C0Mn0N/ABA/BA17 .BLEX 
COMMON/CIJR/noSL (10) 

COMMON/TOL/TOL 1 , T0L2, T0L3 , T0L4 
C0MMUN/TR5/TRA5 

CnnrtUN/GAS/CP >GAM . EMME . ERRE >FXP 1 .EXP2 , UTS2 IS3 

COMMON/CnZ/TCOR, iCnZ, IlNC, TAX , ICL, IAN, ICON 

rOMMDN/MAC/TN 

COMMON/IUIT/IND,INZ,IWR 

OOMMON/AUS/XCL 

COMMON/CSS/C J , G , 01 

CnMMON/UARl/RC( 10) . RS0LD2 , RS0LD3 . RSOL D4 , ASFO . RSFO RRFO.ARFO, 

*RR ( 10 ) >RRnLD2,RRQLf)3,RR0LD4,CU,CR ,UA1 ( 1 0 ) , DALF ( 1 0 ) ,DBET(10 ) > 
)tcASF AMs Bl(20) / ^ / / / 

COMMON/0aR 2/B6(20) .ZR , ZS . ARF > 82 ( 20 ) > PR , AMR , VUl ( 1 0 > 

COMMON/UAR. 1/PTE ( 1 0 5 , RS 1 , RSl , RS5 , T2 ( i 0 ) 

C0MM0N/UAR4/BR1 >BR2:bR3>RR1 , RR3 , RR5 , U A2 ( 1 0 ) 

nOhMON/UARS/PRATK 10) .KINCK 10) ,ALPAli ( 10) ,BETAil( 10) ,ZETAi( 10 ) 

( 1 0 ) , ALFA22 ( i a ) , BET A22 ( 1 0 ) 

C0MH0N/gAR6/PT3( iO ) >TT2(i0) ,PTl(i0) , DELH ( 1 0 ) , ALFA? ( 1 0 ) 'vUJ2(i0) , 
JKWR2< 10) .T2S< 10) ,T?TS( 10) 

COMMON/OAR7/TTIS(5) .BETAT(5) ,ETAT(5) ,ETAI( 10) ,ETAS(10) ,ETAR(10) ; 
*RSTAR(i0) ,AKIS(10) .PSIR(IO) 

CnMMON/UAR8/DRl(10 3 ,AMW1( 10) , AMU2< 1 0 ) >BETET ( 1 0 ) .PR ATIT ( 1 0 ) . 
.1cAhR2( 10) , YS(iO) ,Xl(i0) ,AREAi(10) ,ZETAPS<10) ,WPERi( 10 ) ,YR( lO) , 

n0MM0N/UAR9/ZETAR( 10) .ZETAPRdO) , AS ( 1 0 ) , AR ( 1 0 ) , SI 1 ( 1 0 ) ,SI2< 1 0 ) , 
*S1 ( 1 0 ) , DST)X 1 (10),WT1(10),HE(10) 

CnhhON/UARlO/UUl ( 10) ,DHEDX( 10) , DSDX2 ( 1 0 ) ,R I 1 ( 1 0 ) ,RI2(10) , 

)itRI3( 10) ,RI4( 10) ,RI (1.0 ) ,SR1(1 0) ^SR2(i0 ) 

COMHON/UARll/YOLDdO) » AA d 0 ) . SR ( 10 ) ,PRAT2d0) ,WPER2dO) , 

JKDUJDXdO ) .TllSdO ) ,PPAT3dO ) .88(10) . ALFAdO ) 

COMMON/UAR12/BETA(10) ,D£LR(iO) ,WPERO( 10) 

:*cZETAR3(2n> 

COMMON/UAR^ 



R12/BETA(10) ,D£LR(iO) ,WPERO( 10) .ZETASdO) >ZETAR1 (20) , 
) ,ZETAR5(20) ,R1(20) , Al (20 ) , T1 0 (20 ) 

R 13/ST 1(20) .TRR(20) .R2(20) ,A2(20) , RING (20) ,DRdO), 



JKBET 0 ( 1 0 ) , STAL 1 1 (1 0 ) , f^R EA2 ( 1 0 5 /OR 1 ( 1 6 ) 
COr.MON/UAR 14/ULBM , PR ATS . OMEG 



67 

750 

7750 

100 



530 



COMMON/ALl/ALFAl ( in) ,Ul(iO) ,TTn,PPM.RSd 0) , SI , TNI . H . D ,CI , T1 d 0 ) ; 
XcPl (10 ) .TO.TFI . ALI .BF8P,XX,Ai<^G2n'AMSi (10) ,8 .TN,C. Ti^, AL,S0;TN0, 
>»cCO , TEO , U ( 1 0 ) , 6 i I , D 1 0 , D2I . D2n , ANG2.I . ALE AX , T t > 1 . P t H , ALO . AMC 
COMMON/ AL2/BETA2dO>; BETA! d 0) ,BETA0dn),U2( 10) , TTE (1 0 ) , U2 d 0 ) , 
)«cSIR ,TNIR .HR . DZ , CIR , T IPC , SZ , TNR ,CR ,SQR ,TNCR ,COR . ALT R , At R , ALOR . 
-YP2(iO) ,yn?( 10) ,W1 (10) ,TEIR,TER,TE0R,D1 tR,DiOR,BETAZ,BETAl ,ANM, 
*TIR.TR,TnR,STALI(10) 

COMMON/ARE/REE 

COMMON/ fRA/XPOl (5,8) , XP02 ( 6 , B ) . ALF 1 (8 ) , ALFOl (5 ) , ALF02(6) , 

)^:Yd0) , Yt (10) .Q(6) ,RX(30) ,RY(30) , IR(30 ) ,Z(6) ,C1(4,8) ,C2(4,8) 

DATA NAME /2HP A , 2HR T , 2H2 / 

DATA NAMR/2HPA,2HRT,2H3 / 

DO 67 1=1,5 

U( I )=RPM*5. 14159/30 ./12.*RS( I) 

DR(I)=0, 

DELR(I)=0. 

7ETAS<I)=.iO 
ZETARd ) = . 15 
ZETAPS(I)=0 . 05 
ZETAPR( I)=0 . 05 
COSL(I)=1.0 
YS(I)=1.0 
YR(I)=i . 0 
N9=0. 

N5 = 0 

N9=N9+1 

CONTINUE 

RS(2)=RS0LD2 

R8(3)=R80LD3 

DO 530 I=l>5 

XI (I)=RS(I)/RS<3) 

ASF=ASFO 
RSF=RSFO 
FSi=i . 0 
F82=i . 0 

CALL CHAN ( TTO , AMC , PTO ,RC , ULBM , UCHAN , WPER 0 ) 

NS=0 
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0079 

0030 

0081 

0082 

0083 

0084 

0085 


810 


0086 

0087 

0088 

0089 

0090 

0091 


120 


0 0 92 


801 


0093 

0094 

0095 
0 096 

0097 

0098 

0099 

0100 
0101 
0102 
0103 


802 


0104 

0105 

0106 

0107 

0108 
0109 


71 


0110 

0111 

0112 

0113 

0114 

0115 


201 


0116 

0117 


320 


0118 

0119 

0120 


36 


0121 

0122 

0123 

0124 

0125 


310 


0126 


200 


0127 

0128 
0129 


130 


0130 

0131 


4322 


0132 


5000 


0153 

0134 

0135 


220 


0156 

0157 


221 


0138 

0139 

0140 

0141 

0142 

0143 

0144 


880 


0145 

0146 

0147 

0148 

0149 

0150 


860 


0151 


881 


0152 


822 


0153 

0154 

0155 

0156 


861 


0157 

0158 


894 



DO 80i K = i ,15 
CAi.L STATR ( ALF A1 , XI , TTO . P TO , AH5 , T 1 . P 1 , VI . VA 
•r.DSDXl ,UlJi , PRATl . T ( 15 ' SS , DALF " RSF , DEiLr , CV , , 



VAl . 
2E 



STl ,STP,YS,S1 > 
TAP5,RS'RS{ ,R^3, 



_ iMr3,Ax»wL.i ,r 1 C »rTO,TTFTTG, A5F ,7S,RSF, ASF , 

)*:7R,RSF, r ,U*r.HAN, VAl^UPERi,CODF,iiLBf1,Bi ,ffS,TlPC,ARFA{ ) 



*RS5,ZF1 A5;dR ,2ETA5 . AgSI , NS , VR 1 ) 

CALL AL05i(ZETAS,2ETAPS) 

DO 120 T=l,5 
PTE(I)=PTn 
TTE(l)=TTO 

CALL FLOyR(PRATl . 7ETAPS ,Xi , UI1 ,PTE , 

> r ,U*CHAN. VAl UPER^ 

IF'dN.EQ.l) AfiC=AHC-.01 
IF(IN.EQ.l) GO TO 7750 
TF(CODE-20.) 801.302,801 
CONTINUE 
CONTINUE 
FCl=i. 

FC?=FC1**2 
^!,rF=ARFO 
RRF=RR0LD3 
RR(2)=RR0LD2 
RR(3)=RRF 
RR(4)=RR0LD4 
DO 71 I=i>5 

A-lSl ( I ) = VI ( I ) /SORT ( GAM)*cERRE*G)ltTi ( I ) ) 

X2(I)=RR(.T)/RR(3) 

CONTINUE 

CALL ROTOl (VU1,VA1 .RPH,U. PETAL .HE .TTF, PTE .X2, 

)KRS,ZETAR ,2ETAPR .RR, DHEDX, DSDXl , S 1 , U2 , OHFC . pRi , 

% 7ETAR,PA,RS1 ,RP3,RS5,PET0,STALII,RINCI,VR1) 

CODE=l. 

IhACC=0 
DO 200 K=1 , 14 

CALL ROTO? (BETA2.HE,DHEDX,DSOXi,DSDX2 .VA2.UL12 
JfcYR.ZETAR ,RX1 .RI2.Rl3jRl4,Rt ,SRl,oR2. AA,5R ,TTE, 
*T2 S.INDS:dBET,RRF,DFLR ,CV,CK,DR,RR,RR1 ,P.R3,RR5 
CALL AL0.52(2ETAR,7FTAPR) 

IF (INDS-5 ) 310 .320,310 
WRITE (6,36) (AA( I) , 1=1,5) 

WRITE (6. 36) ( (AA( I ),7FTAR(T ) ,VA2(I) ) ,1=1 ,5) 

FnRMAT(35H ENTROPY INDETERMINATE, PRINT Aa 1-5 
t 1-5,5E12.4/25X,10H VA2 i-5,5Ei2,4/> 

TND=i 

CALL FL0WR(PRAT2.7ETAPR,X2,WX1 . PTE , PTO , TTE . TTO 
*7R,RRF,2,WCHAN,VA2,WPER2,C0DE,UL&rt,B2,RR,TIPC, 

IF(IN.EO.i) AMC=AMC-.01 

IFdN.EQ. 1) GO TO 7750 

CONTINUE 

CONTINUE 

XMACC=IMACC+i 

IFCIMACC.GE. 10) GO TO 220 

FORrtATC /20H LOOP IN SLINE ROT//5E10.3) 

IF(CODE-40.) 201,220,201 

CONTINUE 

DO 221 1=1,5 

DELR (I)=RS(I)-(RC(I)+RR(I))/2. 

DR(I)=RC(T)-RR(I) 

CnSL ( I ) =SORT ( CV*CV/ V DR d ) )|:*2+CV*CV ) ) 

N3=i 

DO 881 K=1 , 15 

CALL ALOFj (ZETAS,2FTAPS) 

CALL STATR ( ALFA< , X , TTO . PTO , AMS , T 1 . P 1 . VI . VAl , ST L . ST 2 , YS , S 1 . ^ 
1JDSDX1 .VU1 , PR ATI , Ti IS SS , 6 aLF ' RSF , DELR , CV , CK , 2ETAPS , R S , RSI , R§3 , 
*RS5.Z£TAS,DR ,ZETA1,AMS1,NS,VR1) 

DO ^60 T=L,5 
PTFd)=PTO 
TTE(I)=TT0 

CALL FLOWR (PPATi,7ETAPS,Xl .WTi . PTF , P TO . TTE , TTO , ASF , 7J 
*ZR,RSF, L ,wrHAN,VAi. ,WPERl,Cf)DE,liJLBM,Bl ,RS,TIPC,AREAi ) 
TF<IN.^g!i) AMd=AMd-.01 
TFdN.EO.i) GO TO 7750 
IF(CODE-20.) 881,822,881 
CONTINUE 
CONTINUE 

CALL ROTOi ( VUl ,VAi ,RPM .U .PFTAl . HE , TTF . PTE , X2 , P 1 , T1 .U1 .UUi ,X1 , 
*RS,7£TAR ,ZETAPR,RR . DHEDX . DSDXi. , S 1 , U2 , OMFO , BR 1 , BR2 , BR3 , FSl , FS2 , 
*7ETAR ,B6,RSi ,RS3,RP5,BET(!),STALII ,RINC,VR1 ) 

CODE=l. 

DO 896 K=L,10 

CALL AL0S2 ( 2ETAR , 7ETAPP ) , 



PI .T1 ,U1 ,UU1 ,X1 , 
BR3,BR3,FS1 ,FS2, 



W2,VU? .V2.X2.U2, 
PTF T2,P2,AraT2, 
,NS,UR2) 



,5E12 .4/,25X,10H2ETAR 
^ASF^2S,RSF,ARF, 



.2S,RSF,ASF, 
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CALL ROTO? (BETA? .HP, .DHPDX ,D5PX1 .PSDX? .UA? .UU? ,UL? .VU2.^? X? .U2, 
-XYR .ZET^R.RTl . I? . R T .1 , fy 1 4 , R T ' SR 1 . SR2 > AA , , TTE > PT“, T? , PRAT^, 

:<T?S,INDS,nBET RRE .DEI R .C(>,CR , DR '.RR ,RRi .RRS ,RR5 ,NS,UR2) 

IF(lNDS-l) 895,320.895' 

895 CALL ELOWR (PRAT2 .ZETAPR ,X2 .UTl . PTE ,PTO , TTE ,TTO . ASF ,7.S,RSF , ARE ,ZR > 
•1.RRF .2,MjCHAN,9A2,WPPR?,C0D£;WLBrt,B2,RR ,TIPC,AREA2) 
lEClN.EQ.l) AhC=AhC-.0l 
IF(IN.EfJ.t) GO TO 7750 
TF<CODF.-?0 . > 896,897,396 
S96 CONTINUE 
897 CONTINUE 

226 DO 227 1=1,5 

227 PRAT3(I>=PT0/P2<I) 

IND=0 

PRAT5=(PPAT3(l)+2 . * ( P R AT3 ( 2 ) +PR AT3 ( 3 ) +PR AT3 ( 4 > > -fPR AT3 ( 5 > ) /8 . 
UlRITEd ,?65> PRAT5 

265 FORhATdX. ^COMPUTED PRESSURE RATI0=SF6.3) 

DIFF=AB.5(PR-PRAT5) 

TAL=10L3*PR 

IF(TAL-DIFF) 920,910,910 
910 Nil=0 

no TO 223 
920 CONTINUE 

710 IF(PR-KRAT5) 712,712,714 
712 AhC=ArtC-DIFF/18 . 
no TO 750 

714 ArtC=AMC+DIFF/18.- 
no TO 750 
223 CONTINUE 

CALL £XEC(8,NAMR) 

END 
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PROGRAM PART3(5) 

DIMENSION NAME<3) 

DIMENSION NAHR<3) 

CfJhnaN/ARA/BA17,BLEX 
COMMON/CUR/nOSL( 1 0 ) 

GOMnON/ THI /TOL 1 , TOI .2 , T0L3 , TOL 4 
C0MM0N/TR5/TRA5 

COMMQN/GAS/CP ,GAM >FMME.ERRE.EXP1 ,EXP2,VTS2,VIS3 

COmMON/CDZ/ICOR , icnz, IlNC, iAi , ICli , ian , icon 

nOMMON/MAC/IN 

CnMMON/IUI/IND,INZ, IWR 

COMMON/AUS/XCL 

COnMON/CSS/CJ,G.01 

CnMMON/yARl/RC( iO ) .RSnLD2.RSnLD3.RS0l D4 , ASFQ . RSFO . RRFO > ARFO , 
J^tRR(lO) >RRC)LD2.RR0LD3,RR0LD4,Cy,CK,MAl(1.n) , DaCf ( 1 0 ) , DBET ( 1 0 ) , 
:KASF..AMd,Bl(20) 

COMMON/yAR2/B6(20) ,7R . ZS . APF , B2 ( 20 ) .PR , AMR ,UU1(10) 
C0MHGN/'JAR3/PTE ( 1 0 5 . RS3 ,I^S5 , T2 ( i 0 ) 

COMMON/UAR4/BR1 >BR2 .BR3.RR1 ,RR3,RR5.GA2(10) 

C0MM0N/GAR5/PRA1 1 Cl A ) .RlNCIClO) ,ALF AlK 10) , BET A1 1 ( 1 0 ) ,ZETA1 ( 10 ) , 
.*V2(10) ,ALFA22(10) .BFTA22(10) 

CCMMON/UAR6/PT2(l0) ,TT2(10) , P T1 ( 1 0 ) , DELH ( 1 0 ) , ALFA2 ( 1 0 ) , GU2 ( 1 0 ) , 
*UIR2(10) ,T2S(10 ),T2.TS(10) 

CaMMQN/'7AR7/TTISI5) ,BETAT(5) ,ETAT(5) ,ETAI ( 10) ,ETAS( 10) .ETAR ( 10) > 
*RSTAR(i0) .AKIS(IO) .FSTRCIO) 

C0hMQN/UAR8/DRl<in) .AMUKIO) .AMV?(10) >BFTFT(10) .PRATIT(IO) . 
)KAMR2(10> >YS(10) ,Xl(iO),ARF.A1. (10) ,ZETAPS(10) ,WPER1(10),YR(10) . 
«X2(10) 

C0MM0N/GAR9/ZETAR( 10) .ZETAPR(IO) , AS< 10 ) ,AR ( 10 > ,SI1 ( 10 ) ,SI2( 10 ) , 
JkSKlO) ,DSDX1(10) ,WIl(i0).HE(10) 
rOMMON/GAPiO/UUl (10) DHFDXO 0) , DSDX2 ( 1 0 ) , R T 1 ( 1 0 ) ,RT2(i0) , 
:icRI3(l0) .RT4(10) .RI(1(^) .SRl(lO) ,SR2(10) 
nOMMON/OARll/YOLD(i 0 > . AA( 10 ) ,SR C 10 ) .PRAT2( 10) ,UPER2(10) , 
JtcDWDXdO) .TIIS(IO) .PRAt3(10) , SS ( 1. 0 ) . Al FA ( 1 0 ) 
r:QMMON/UAR12/BeTA(l(i) ,DELR(1 0) , WPE R 0 ( 1 0 ) , ZETAS < 1 0 ) , ZETAR 1 < 20 ) , 
)K7ETAR3(20) ,ZETAR5(20) f R l ( 20 ) , Al ( 20 ) , T 1 0 ( ^ 0 ) 

GOMMON/VAP.i3/STl(20) ,tRR(20) ,R?(20) .A2(20) .RING (20) ,DR(10) , 
*PFTO( 10 ) .STALIKIO) ,AREA2( 10) ,UR1( 10) 
CnrtM0N/UAR14/ULBf1,PRAT5,0MEG 

COMMON/ALl/ALFAK 10) ,‘J1(10) .TTO.RPM.RSdO ) , SI, TNT .H.D,CI,T1(10), 
.'KPldO ) ,TO,TEI,ALI,BFSP ,XX, ANG20;AMSid 0) ,S'TN,C,TE. AL,SO'tNO, 
*CO,TEO,U(iO) ,D1I .DiO,D2I .D20,ANG21 .ALFAX.TT ,T .PTO , ALO , AMC 
COMMON/ AL2/BETA2dO) .BETAldO) .BETAOd 0) . W2 (1 0 ) , TTE (1 0 ) , 02 d 0 ) > 
*SIR,TNIR .HR .DZ,CIR,TlPC,SZ,TNR ,CR,S09 , TNOR . COR , ALI R , ALR ,ALOR . 

, W02 ( i 0 ) , W1 ( ( 0 ) , TEIR , TER , TEOR , Di TR , D1 OR , BETAZ , BETAl , ANM > 
*TIR»TR,TOR,STALldO> 

COMMON/ARE/REE 

C0MM0N/TR.A/XPQ1(5,B) , XP02 ( 6 , 8 ) , ALF 1 ( 8 ) , ALFOl ( 5 ) , ALF02 ( 6 ) , 

)|cY(10) . YldO) ,0(6) ,RX(30) ,RY(30) ,IR(30) ,Z(6) ,C1 (4,8) ,C2(4,a> 

DATA Name /2HPA,2HRT,2H2 / 

DATA NAMR/2HPA,2HRT,2H3 / 



PAGE 



RPM 



TOTAL/STATIC 



PRESSURE RATIO 



999 FORMATdHl) 

URITE(6,999) 

URITE<6,401> 

401 F0RMAT(///27X, ' SET 
XTOTAL INI ET TOTAL d 

UPITE(6,402) 

402 FORMAK 27X . 67HNUMBER NUMBER 

)KRE TEMPERATURE) 

403 F0RMAT(72X, 5H(PSI), 7X,8H(DEG, R)/) 

J=1 

IS=1 

U1RITF(6,4 05) J,IS,RPM,PR .PTO.TTO 

405 FORMAT (27X, 13, 18, FU . 1 ,F14 , .3,F14 , 3,F15 . 2/) 

UiRlTE(6,404) 

404 FORMAT (//57X21H STATOR EXIT SOLUTION///) 

WRtTE(6,406) 

406 FORMAT(iX. ^STREAM RADIAL X=K/RM RADIAL BLADE 

*/UArt BLADE LOSS ZETA* FLOW RATEd 

URITE(6,4D7) 

407 FORMAT(iX>109H LINE POSITION SHIFT OPENING 

EFFICIENCY COEFFICIENT CONTINliTTY FRACTION /) 
WRITE (6, 411 ) 

411 F0RMAT(12X,4H(IN) , 13X , 4H ( IN ) , 5X , 4H ( IN ) ) 
on 408 1^1,5 



INLET 



PRESSU 



Y=VA 
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0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

on SB 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 



PPATi <T)=1 ,/PRATl(I) 

. ALFAil (T)^AI..FAl(I))<c57.3 
SSTAtl ( I )=r<FTAl ( I )*57 . 3 
ALFA2C I'^=ATAN( VU2( I )/VA2< I ) ) 

M2( 1)=VA?( I)/C0S(ALFA2(I) ) 

U2( 1>=BQRT(V2(1 >¥V2( I )+US2( I ))KWR2(I) ) 

ALF A22 < I ) =ALr A2 ( I ) )K57 . 3 
BETA22( I)=BFTA2( I))K57 .3 

DFLH < I ) = ( U < I ) Wl\ 1 ( I ) **112 ( I > *VU2 < I ) ) / ( G^C J ) 

TT2 < I ) =TTO-DFLH < T ) /CP 
PT2(I)=P?(X))tc( ]T2( I)/T2a ) )*){(eXPl 
PTia)=Pl(I)^< TTO/Tl ( T ) ) r^EXP 1 
T?IS< J )=TT0Xi(H2< I )/PT0))|cicF.XP2 
Tr’S( I )=TTF( I)>K(P2( I)/PTE(.T) ).k*eXP 2 
TTISd )=TTn)»i(PT2( I>/PT0>)|r^EXP2 
BFTAK I)=PT0/PT2< I) 

£TAT(I)=(TT0-TT2( I) )/(TT0-TTT8(I)) 

ETAKI ) = (TT0-TT2(I) )/(TT0-T2TS(I) ) 

ETAB(I) = (TTn-i id) )/(TTO-TlJB(I) ) 

FTAR d ) = <TTE( I)-T2( T) )/(TTE(I )-T2S( I) ) 

P3TAR(I) = (TlIBd )-T2TSd) )/<TT0-T2ISd) ) 

AK T S ( I ) =2 . J*CP ( TT0-T2IS d ) ) /IJ d ) 

PBJRd)=BqRT(b'.TAR(I)) 

DRi < I )=RC(I )-RS< I) 

AMWld)=Wl(I)/SQRT(GAh)HERRF*r,)lcTid) ) 

Ai'iBi ( I ) =yi ( T )/SQ:?T(GAii^ERREiKG)lcTi ( I ) ) 

AMV2d)=^2d )/SQRT(GArt*EkRF)*cG*T2(I) > 

AMR2 ( I ) =U2d ) /SORT < RAh)*cERRE*G*T^?( I ) ) 

BETET(I)=PTE(1)/P2d) 

PRAT1T(I)=PT0/PT1 ( I) 

40 8 WRITE (6 .409) 1 , k S ( I ) , XI (I ) , DR 1 < I ) , AREAl ( I ) ,YSd) , ETAS ( I ) , ZETAi ( I ) , 
)lc7.F.TAPS(I),WPERid) 

DFL.H(iO) = 0. 

DO 240 T=i,4 
L=I + 1 

240 DELH(i0)=DELH(10)+.5:*c(WPER2(L)-UPER2d) )5 *c(DELH(L)+DELH(I) ) 

HP=DELH< 10 )*CJ*WLBM/550 . 

AhOrt=HP*5BO ./OHEG 
THETA=SQRT(TT0/518. 4) 

DELTA=PT0/14 .7 
HP1=HP/(THFTA*DELTA) 

AMOhi=ArtOH/DELTA . 

RPMi=RP«/THETA 
UI.BMi=WLBn>|{THET A/DELTA 

ET A5= < ETAI ( 1 ) -t-ET A I < 5 ) +2 . * ( F.T AI < 2 ) +ET AT ( 3 ) 4-RTAT < 4 ) ) ) /8 . 
BETA6==(BFTAT( 1 )+BFTAT(5)+2 . :ic< SETAT(2 >+PFTAT(3)+BFTAT( 4) ) )/8. 
ETA6=(ETAT(i )^ETAT<5)+2.*(FTAT(2)*»-FTAT(3)+ETAT(4) ) )/8 . 

AKTS5=( AKIB( 1 )-4-AK IS < 5 ) +2 . 3*^ < AK IB ( 2 ) +AK TB ( 3 ) +AK IS ( 4 ) ) )/8. 
RSTAR5=(RBTAR(i )+RSTAR(5)+2.:*c(PBTAR (2)+RSTAR(3)+RSTAR(4) ) )/8. 

409 FnRMAT(iXd4,Fi2.3,Fi0.3,F9.4,F9.4,Fli .4,Fii,4,Fl4. 4,Fi3,4,F14.4) 
WRITE(6>4i2) 

412 F0RMAT(///22X,23HABS0LUTE VELOCITY ( FPS) ,27X ,23HRELATIVE VELOCITY 
3|c(FPS)//) 

413 FORhAT(iX,6H3TREArt.2X,2(50H AXIAL RADIAL TANGENTIAL OVER 

MLL ) ,7H WHEEL ) 

URITE(6,413) 

WRIT£(6.414) 

414 F0RHAT(iX,6H LINE , 2( 50HCOMPONENT COMPONENT COMPONENT VE 

3»:L0CITY ) ,8HVELnCITY//) 

DO 415 1=1,5 . 

415 WRITE<6,416) I,VA1(I) ,VRld) ,VU1(1) , VI Cl ) , VAl ( I ) , VR 1 d ) ,WU1(I) ,U1( 
* I ) , U d ) 

!,3F12.2,6X) ,F8.2) 



MACH NUMBER 

PRESSURE 



(PSD 



FLOW ANGLE 

PRESSURE) 



(DEG) 



416 F0I^MAT(I5,2X,2(F8. 

WRITE(6,4i8) 

418 F0RMAT(///7X,il3H 

* TEMPERATURE 
URITE(6,419) 

419 FORhAT( 7X.113H 

* (DEG. R) 

WRITE(6,492) w 

492 FQRMAT( 7H STREAM) 

WRITE(6>420) 

420 FORMAT (7H LINE ,2(24H ABSOLUTE RELATIVE ) ,2(24H 

* STATIC ),24H TOT/TOT TOT/STA /) 

DO 422 1=1 >5 

422 WRITF(6.42i )I,AMSi (D ,AMU1( I ) , ALFA! 1 ( I ) , BETAil (I ) , TTO ,T 1 ( I ) , PTl ( I ) 
*,P1 (I ) .pRATlTd) .PRATi d) 

421 FORMAT(I4 ,3X,2Fin .?,4X,2F10 .2,4X,2F10 .2,4X,2F10 . 3 , 4X , F 1 1 . 4 , FI 0 . 4> 
WRITE (6, 9^9) 

WR lib (6,401) 



RATIO /) 



TOTAL 
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0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 
0193 

0184 

0185 
0136 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 
0201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 
0227 
0229 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 



ROTOR EXIT SOLUTION///) 



URTTE(6,402> 

■ URITE(6.403) 

452 F0RMAT(//57X,21H 
IS=2 

URITE( 6,405) J , 1 S , RP h , PR , PTO , TTO 
UiRIT£(6,452) > ^ ^ 

WRT.TE(6,406) 

UlftTTE(6;:407) 

00 423 1=1,5 

WRITE(6,4i2) 

WRITE(6,413) 

WRITE(6,414) 

DO 424 1=1,5 

424^WRITE(6,4j6)I,UA2( I) ,WR2< I) ,UU2(I) ,U2(I)*,UA2(1) ,WR2(I) ,UU2(I) ,U2(I 

Wl^ltE(6,418) 

WRITE(6',419) 

WRITE:<6,492) 

URTTE<6' 420) 

DO 425 1=1,5 

425 WRlTt(6,42l) I ,AMU?(I) , AMR2 ( I ) , ALF A22 ( I ) , BETA22 ( I ) ,TT2( I ) , T2( I) , 

*PT2(I),P2<I) , BETAT(I) ,PRAT3(I) > > > 

URITE(6,491) 

491 FORMAT(///) 

WRITE(6,426) 

426 FORMAT (7H STREAM, 41H 



EQUIVALENT 

TEMPERATURE 



EQUIVALENT 

INLET 

PRESSURE 

(PSD 



427 FnRMAT(7H LINE , 38H 

428 FORMATC 7X, 38H 

429 FORMATC 7X , 22H (DEC;. R) 

URITE(6,427) 

WRITE(6,428) 

WRITE(6,429) 

00 430 1=1,5 

430 WRITE(6,431 )I,TTE(I) >PTE(I) ,BETET(I) 

431 F0RMAT(T4,Fi3.2,F15.3,Fli .1) 

WRITE<6,999) 

WRITE(6,401) 

WRITE(6,402) 

URITE(6,403) 

IS=3 

14RITE(6,405) J , IS , RPM , PR ,PTO ,TTO 
UIRITE(6,441) 

441 FnRMAT(//45X,31H0VERALL TURBINE CHARACTERISTICS///) 
URITE<6,442) 

“ STREAM PRESSURE RATIO 

BLADE/ JET THEORETICAL 



EQUIV/STATIC) 

PRESSURE) 

RATIO) 



LINE 

SPEED 



) 



TflT/STA TOT/TOT TOT/STA 
RATIO DEGREE OF REACTION /) 



TOT/TOT 



442 FORMATC i02H STREAM PRESSURE RATIO EFFICIENCY 

^ head ■ 

WRITE(6,443) 

443 FORMAT <in2H 
^ COEFFICIENT 

DO 444 I=i ,5 
BLAJF=1 ./SQRT(AKIS(T) ) 

444 WRITE (6, 445) I , PR AT3 ( I ) ,BETAT( I ) , ETAI C I ) , ETAT ( I ) , AK IS C I ) , BLAJE , RST 

^AR (I)^ ' ^ ' 

445 F0RrtAT(I5,F14.4,Fll.4,Fii .4 ,F13.4,F12 .4,F15.4,F16,4) 

UIRITE(6,446) 

446 F0RMAT(///53X.24HMASS AVERAGED QUANTITIES///) 

447 FnRMATC52X,l3l^l-:0RSE POWER =, FI 0 . 2 , 3X, 4H(NP ) ) 

448 FORMAT (52X, 13HM0MENT = , F I 0 . 2,3X , 7H ( FT-LB ) ) 

449 FnRMAT(S2X, l3HFLnU RATE =, FlO . 2 ,3X,BHU.B/SEC)//) 

461 F0RMAT(43X,22HREFFRRED RPM =,Fi0,2) -v 

462 F0Rr1AT(43X,22HREFFRRFD HORSE POWER =, F 1 0 . 2 , 3X , 4H C HP ) ) ‘ 

463 FORMAT (43X,22HRFFFRRF.D MOMENT = , F 1 0 , 2 , 3X ; 7H ( FT-LB ) ) 

464 F0RMAT(43X:22HREFFRRE0 FLOW RATE = , FI 0 . 2 , 3X , 8H < LB/SEC ) // ) 

465 FnRMAT(4nx;25HT0TAl., /STATIC EFFTCIENCY =,F10.4) 

466 FnRMATC40X,25HTOTAI. /TOTAL EFFTCTENCY= , FlO. 4) 

467 FORMAT (36X,;?9HT0TAI. /STATIC PRESSURE RATtO =,F10,4) 

468 FnRflAT(36X.29HTDTAL/T0TAL PRESSURE 

469 FORMAT (34X;31HHFAD COEFFICIENT 
471. FnRMAT(34X^3iHTHE0RF7TCAL DECREE OF 

472 FORMAT (34X , 3) HBLADF/JFT SPFF.D RATIO 

473 F0RMAT(34X .31HMACH NUMBER AT STATION 
WRITE(6,447) HP 
WRITE(6,448) AMOM 

yj 

RPMl 
HPl 
AMOMi 
WLHH1 



RATtO =,Fin .4//) 
=,F10 .4) 
REACTION=,F10 . 4) 
=,Fin .4) 
0 =,F10.4) 



WRITt(6,449) 
WRITE(6,46l ) 
WRITfc (6,462) 
LR ITEC6;463) 
W?TTFC6.464) 
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0.?39 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 
024R 

0249 

0250 

0251 



WRTTE<6, 
' WRT1E<6, 
5.mRIT£(6, 
WRITE (6, 
IaIR I TE < 6 , 
MLAJFb=i 
MRITE(6, 
WRITt(6, 

wri:te(6, 

IF(INZ-1 
930 TF(NU-i 
400 CONTINUE 
END 



465) ETA5 

466) ETA6 

467) PRATS 
4cv8) BETA6 
469) AKIS5 
./5QRT(AKIS5) 

472) BLAJE5 
471) RSTAR5 

473) AhC 

) 400,930,930 
) 400,400,400 
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